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The  influence  of  surfactants  on  the  motion  of  bubbles  and  drops  in  a  Hele-Shaw 
cell  is  studied  by  carrying  out  experiments  in  a  variety  of  fluid-fluid  systems.  Saffman- 
Taylor  theory  predicts  the  translational  velocities  and  shapes  of  drops  and  bubbles,  not 
influenced  by  surfactants.  In  order  to  differentiate  the  cases  with  and  without  surfactants, 
the  motion  of  air  bubbles  and  water  drops  driven  by  silicone  oil  was  studied.  A 
controlled  amount  of  surfactant  is  added  to  the  water  drops  in  order  to  study  its  influence 
on  the  bubble  motion.  In  the  case  of  air  bubbles  in  silicone  oil,  surfactants  were  likely  to 
have  negligible  influence,  although  they  may  be  present  as  contaminants.  As  expected,  air 
bubbles  were  moving  much  faster  and  their  shapes  were  also  in  accordance  with  the 
available  theories  for  the  surfactant  free  system.  On  the  other  hand,  surfactant  containing 
water  drops  moved  with  much  smaller  velocities  and  their  shapes  were  also  distorted. 
These  observations  are  apparently  consistent  with  those  of  Kopf-Sill  and  Homsy.  When  a 
surfactant  containing  drop  is  set  in  motion,  the  resulting  spatial  variation  of  surfactants  at 


the  interface  induces  Marangoni  stresses  along  the  interface,  and  consequently  the  drop 
surface  may  become  rigid.  If  the  surrounding  fluid  is  the  wetting  phase,  the  net  drag  force 
experienced  by  the  drop  in  the  wetting  film  region  is  significantly  high  resulting  in 
smaller  drop  velocities.  Because  the  surface  tension  has  a  spatial  variation,  to  keep  the 
pressure  inside  the  drop  constant,  the  curvature  adjusts  itself  resulting  in  a  variety  of 
bubble  shapes.  Assuming  that  the  bubble  surface  is  rigid  due  to  surfactant  influence,  the 
translational  velocity  of  an  elliptic  bubble  is  estimated.  The  result  indicates  that  the 
bubble  velocity  can  decrease  by  an  order  of  magnitude  compared  to  the  prediction  of 
Taylor  and  Saffman.  The  retarded  bubble  velocity  is  apparently  in  reasonable  agreement 
with  the  experimental  observations  suggesting  that  the  puzzling  observations  of  Kopf-Sill 
and  Homsy  are  likely  to  be  due  to  surface  active  contaminants. 

Finally,  two  phase  displacement  of  a  shear  thinning  power-law  fluid  by  an 
inviscid  fluid  is  analyzed.  The  flow  domain  is  divided  into  three  different  regions  and  a 
leading  order  solution  to  the  flow  field  and  the  interface  shape  is  obtained  using  the 
technique  of  matched  asymptotic  expansions.  In  the  case  of  Newtonian  fluids,  capillary 
number  uniquely  determines  the  film  thickness.  On  the  other  hand,  in  the  case  of  shear 
thinning  power-law  fluids,  the  relationship  is  unique  to  a  given  fluid  and  flow  conditions. 
Finally,  the  flow  field  away  from  the  interface  in  the  power-law  fluid  region  is  also 
solved  assuming  that  the  gradients  of  components  of  velocities  is  negligible  in  the 
directions  normal  to  the  gap  direction.  Unlike  the  Newtonian  case,  the  depth-averaged 
velocity  field  is  found  be  non-Darcian.  A  stability  analysis  is  also  carried  out  to  determine 
the  wave  number  with  the  maximum  growth  rate.  The  result  indicates  that  the  interface 
will  be  less  sinuous  for  the  power-law  fluid  than  the  Newtonian  fluid. 
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CHAPTER  1 
INTRODUCTION 

1.1  Background 

Hele-Shaw  cell  has  been  used  as  a  model  system  to  study  various  physical 
phenomena  involving  multiphase  flows  in  porous  media.  Although  the  concept  of  Hele- 
Shaw  cell  is  more  than  a  century  old,  there  has  been  a  lot  of  interest  in  recent  years  in  the 
motion  of  two  phase  flows  in  a  Hele-Shaw  geometry.1"15  The  displacement  of  a  more 
viscous  liquid  using  a  less  viscous  liquid  and  the  related  fingering  phenomena  have  been 
extensively  studied  both  experimentally  and  theoretically,  using  this  model.2,14"15  When 
the  fluid-fluid  interface  is  a  closed  contour,  the  inner  fluid  is  referred  to  as  a  drop  or 
bubble  depending  upon  whether  the  enclosed  fluid  is  a  liquid  or  gas. 

A  typical  Hele-Shaw  cell  consists  of  a  thin  gap  between  two  flat  plates.  The  gap 
thickness,  which  is  of  the  order  of  a  millimeter,  is  very  small  compared  to  the  dimensions 
of  the  plates.  The  gap  is  filled  with  a  liquid  into  which  a  secondary  fluid  is  injected  to 
form  a  drop  or  bubble  and  its  motion  can  be  studied.  The  injected  drop  or  bubble  is 
squeezed  between  the  two  plates  and  assumes  a  circular  planform  when  it  is  stagnant 
because  of  the  interfacial  tension.  If  the  surrounding  fluid  is  the  wetting  phase,  a  thin  film 
exists  between  the  bubble  surface  and  the  solid  boundary  of  the  Hele-Shaw  cell.  The 
bubble  can  be  set  in  motion  either  by  applying  a  pressure  gradient  in  a  horizontal  Hele- 
Shaw  cell  or  by  imposing  a  buoyancy  force  by  tilting  the  cell.  If  the  characteristic  length 


(a) 


Wetting  film 


7 — 

10  film  » 


5S  (h;gap 


Stagnation  Points 


(b) 


Figure  1.1       Schematic  of  a  moving  bubble  in  a  Hele-Shaw  cell  :  (a)  View  in  the  xy 
plane,  (b)  Bubble  in  the  xz  plane 
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scale  of  the  bubble  in  the  transverse  direction  is  larger  than  the  gap  thickness  (a/h  »1, 
Figure  1.1),  the  dynamics  of  the  bubble  motion  can  be  described  by  Hele-Shaw 
equations16  which  are  a  two  dimensional  description  of  the  gap-averaged  velocity  and 
pressure  fields.  This  approximation  is  valid  over  the  entire  domain  expect  for  a  thin 
region  (of  the  order  of  gap  thickness)  surrounding  the  rim  of  the  bubble.  The  depth 
averaged  equations  resemble  the  Darcy  approximation  of  the  governing  equations  for  the 
flow  through  porous  media.  Hence  the  Hele-Shaw  cell  can  also  be  used  to  study  flow 
through  porous  media. 

Taylor  and  Saffman1  have  solved  the  Hele-Shaw  equations  for  the  motion  of 
bubbles  in  a  thin  gap  with  no  side  walls  assuming  negligible  surface  tension.  Their 
solution  predicts  the  bubble  shape  as  a  function  of  the  area  and  bubble  velocity. 
Furthermore  they  have  concluded  that  for  a  freely  suspended  bubble,  the  bubble  velocity 
(U)  is  always  greater  than  the  surrounding  fluid  average  velocity  (V)  and  the  steady 
shapes  are  ellipses.  Also  if  U  <  2V  these  elliptic  bubbles  travel  with  their  minor  axis 
parallel  to  the  flow  direction  and  if  U  >  2V  they  travel  with  the  major  axis  parallel  to  the 
flow  direction.  If  U  =  2V  the  bubble  shapes  are  circles.  Tanveer3  had  incorporated  the 
effect  of  finite  surface  tension  to  resolve  the  degeneracy  of  the  Taylor-Saffman  solution 
and  obtained  three  solution  branches.  A  stability  analysis4  of  these  three  solution  branches 
indicate  that  only  one  of  them  (McLean-Saffman  branch)  is  stable.  The  inclusion  of 
interfacial  tension  in  the  theoretical  analysis  does  not  alter  Taylor-Saffman's  conclusion 
about  the  relative  velocities  of  the  bubbles  and  the  surrounding  fluid. 


Bubble  motion  was  experimentally  studied  be  several  researchers  using  a  variety 
of  fluid  combinations.5"8  Maxworthy5  had  experimentally  studied  the  motion  of  air 
bubbles  in  an  oil  filled  Hele-Shaw  cell  under  the  influence  of  gravity.  At  smaller  values 
of  the  inclination  angles,  stable  circular  bubbles  or  elliptic  bubbles  elongated  in  the  flow 
direction  are  observed.  However  at  higher  angles  of  inclination,  bubbles  became  unstable 
with  severe  shape  distortions.  When  the  bubble  size  relative  to  the  transverse  cell 
dimension  (a/L)  was  larger  than  0.1,  the  bubble  velocity  was  slightly  larger  than  the 
prediction  of  Taylor  and  Saffman.'  With  decreasing  bubble  size,  however,  the  bubble 
velocity  decreased  to  a  value  smaller  than  the  predicted  one.  This  discrepancy  in  the 
bubble  velocity  was  apparently  resolved  by  including  the  viscous  dissipation  around  the 
bubble  edge  in  calculating  the  bubble  velocity. 

1.2  Motivation 

In  1988,  Kopf-Sill  and  Homsy6  conducted  similar  experiments  to  investigate  the 
shapes  and  velocities  of  various  size  bubbles  driven  by  a  pressure  gradient.  They 
investigated  the  motion  of  air  bubbles  in  a  horizontal  Hele-Shaw  cell  filled  with  a 
glycerin-water  mixture.  Unlike  the  observations  of  Maxworthy,  a  variety  of  unusual 
bubble  shapes  were  observed  as  shown  in  Figure  1.  2.  Besides  the  near-circular  and 
elliptic  bubbles,  "Tanveer"  bubbles,  ovoids,  long-  and  short-tail  bubbles  were  observed. 
By  tracing  the  locations  of  the  leading  and  trailing  edges  of  the  bubbles  for  an  extended 
period  of  time,  they  showed  that  these  unusual  shapes  were  not  transient  but  steady. 
Furthermore,  they  found  that  at  low  capillary  numbers  the  translational  velocities  of  the 
bubbles  (near-circle,  flattened  or  elongated  ellipses,  Tanveer)  were  much  smaller  than 
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Figure  1.  2  Various  Bubble  shapes  observed  by  Kopf-Sill  and  Homsy:  (a)  Near-Circle; 
(b)  Flattened  bubble;  (c)  Tanveer  Bubble;  (d)  Ovoid;  (e)  Short-tail  bubble; 
and  (f)  Long-tail  bubble 


expected.  Especially,  the  velocity  of  near-circular  bubbles  was  about  an  order  of 
magnitude  smaller  than  the  prediction  of  Taylor  and  Saffman. 

In  Figure  1.3,  the  velocity  data  taken  from  Kopf-Sill  and  Homsy5  are  given  for 
three  different  bubble  sizes.  The  capillary  number  in  the  x-axis  is  defined  as  p.V/a  where 
u  and  V  are  the  shear  viscosity  and  the  average  velocity  of  the  glycerin-water  mixture, 
respectively,  and  a  is  the  surface  tension.  When  Ca  is  small,  interfacial  tension  is 
dominant  resulting  in  a  circular  bubble  shape.  According  to  the  theoretical  prediction,' 
the  velocity  of  a  circular  bubble  relative  to  the  average  fluid  velocity  (i.e.,  U/V)  should  be 
2.  Contrary  to  the  prediction,  the  relative  bubble  velocity  was  observed  to  be  only  about 
0.2.  The  very  low  bubble  velocity  at  a  small  capillary  number  was  common  for  all 
bubbles  regardless  of  their  size.  As  Ca  was  increased,  the  relative  bubble  velocity 
increased  slightly  until  it  jumped  suddenly  to  a  much  higher  value  at  a  certain  Ca.  The 
sudden  jump  in  velocity  was  apparently  accompanied  by  a  drastic  change  in  bubble 
shape.  While  the  similar  trend  was  observed  with  all  bubbles  of  different  sizes,  the 
sudden  transition  in  the  bubble  shape  and  velocity  occurred  at  a  lower  value  of  Ca  for 
larger  bubbles.  They  also  reported  a  hysteresis  in  the  transitional  behavior. 

It  should  be  noted  that  the  Taylor-Saffman  solution  exists  only  when  the  relative 
bubble  velocity  U/V  is  greater  than  1 .  Although  the  possibility  of  U/V  <  1  cannot  be 
ruled  out  for  non-zero  surface  tension,  Saffman  and  Tanveer'2  point  out  that  such  a  low 
velocity  of  small  bubbles  cannot  be  explained  by  a  simple  inclusion  of  the  surface  tension 
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effect  and/or  the  wetting  film  effect.  They  suggested  that  the  low  bubble  velocity  might 
be  due  to  a  moving  contact  line  dragging  along  the  cell  plates.  Using  an  ad  hoc  boundary 
condition,  they  showed  that  a  moving  contact  line  could  account  for  the  very  low  bubble 
velocity.  In  Kopf-Sill  and  Homsy's  experiment,  however,  the  liquid  phase  wetted  the 
solid  surface  and  wetting  films  were  always  observed.  Therefore,  the  presence  of  a 
moving  contact  line  was  unlikely. 

We  may  note  that  the  surprising  variety  of  the  steady  bubble  shapes  and  the 
striking  discrepancy  in  bubble  velocities  were  observed  only  in  a  pressure  driven  flow. 
Unlike  Kopf-Sill  and  Homsy,  Maxworthy  did  not  observe  any  unexpected  shapes  but 
only  elliptic  bubbles  elongated  in  the  flow  direction.  In  his  experiments,  the  bubble 
velocities  were  also  in  disagreement  with  the  Taylor-Saffman  theory.  The  disagreement, 
however,  was  not  as  great  as  in  the  Kopf-Sill-Homsy  experiments  and  it  could  be 
explained  by  the  inclusion  of  viscous  dissipation  at  the  bubble  rim.  In  the  work  of  Kopf- 
Sill  and  Homsy,  the  experiment  of  Eck  and  Siekmann7  is  quoted  in  which  the 
translational  velocity  of  small  bubbles  were  also  much  smaller  than  the  predicted  value. 
These  observations  apparently  suggest  that  the  pressure-driven  and  buoyancy-driven 
flows  may  be  intrinsically  different  although  Saffman  and  Taylor2  pointed  out  the 
equivalency  between  the  two  flows. 

The  experiments  of  Kopf-Sill  and  Homsy  involve  several  observations  that  are  in 
need  of  explanation.  These  observations  include  (1)  the  velocity  of  small  bubbles  which 
appear  to  be  an  order  of  magnitude  smaller  than  the  theoretical  prediction,  (2)  a  surprising 
variety  of  steady  bubble  shapes,  and  (3)  the  sudden  transition  of  bubble  velocity  and 


shape  at  a  certain  capillary  number.  Johann  and  Siekmann8  have  measured  the  velocities 
of  air  bubbles  in  a  water  system  containing  steric  acid.  Their  reported  velocities  also 
appear  to  be  significantly  retarded  compared  to  the  theoretical  prediction.  In  all 
experimental  studies  described  above,  the  Hele-Shaw  cell  was  supposed  to  be  clean  and 
free  from  any  contaminants.  Many  fluid  flow  systems,  however,  are  susceptible  to 
contamination  and  it  seems  that  the  perplexing  observations  of  Kopf-Sill  and  Homsy  may 
be  explicable  if  we  adopt  the  notion  of  the  influence  of  surface-active  contaminants  on 
bubble  motion.  One  of  the  objectives  of  this  study  is  to  provide  a  plausible  explanation 
for  the  observations  of  Kopf-Sill  and  Homsy  in  the  context  of  surfactant  influence. 

1.3  Dissertation  Outline 

In  order  to  investigate  the  influence  of  surfactants  on  the  bubble  motion,  we 
conducted  experiments  using  air  bubbles,  water  drops  and  fluorocarbon  oil  drops  in  a 
Hele-Shaw  cell  filled  with  silicone  oil.  Chapter  2  presents  our  experimental  findings  and 
compares  them  with  those  of  Kopf-Sill  and  Homsy.  Air  bubbles  were  chosen  since  the 
surface  properties  of  an  air-oil  interface  were  less  likely  to  be  affected  by  surfactants 
which  might  be  present  as  contaminants.  It  is  well  known  that  unlike  the  surface  tension 
of  water  or  the  interfacial  tension  between  water  and  a  hydrocarbon,  the  surface  tension 
of  non-aqueous  solutions  is  not  affected  significantly  by  surfactants.17  In  case  of  water 
drops,  a  predetermined  amount  of  surfactant  (sodium  dodecyl  sulfate)  was  added  so  that 
the  surfactant  influence  could  be  studied  systematically.  The  motion  of  the  water  drops 
were  found  to  be  drastically  retarded  by  the  surfactant,  whereas  the  air  bubbles  were 
moving  much  faster  as  we  anticipated. 
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When  surface  active  substances  are  present  in  the  system,  they  adsorb  near  the 
fluid-fluid  interface.  Under  dynamic  conditions,  depending  on  the  surfactant  exchange 
kinetics,  a  spatial  variation  of  surface  concentration  can  arise  resulting  in  Marangoni 
stresses.18  As  a  consequence,  the  interface  can  behave  as  rigid  surface.19,20  Hence  in  the 
thin  film  region,  the  drag  increases  significantly  causing  severe  retardation  of  the  bubble. 
Chapter  3  presents  a  theoretical  calculation  that  predicts  the  extent  of  retardation  and 
provides  an  analytic  expression  for  the  drop  velocity  at  low  surfactant  concentrations.  As 
the  bulk  surfactant  concentration  is  increased,  Marangoni  effects  are  expected  to  be 
weaker  and  a  complete  remobilization  is  anticipated.21  Contrary  to  this  expectation,  the 
motion  of  a  drop  containing  a  high  concentration  of  surfactant  is  also  found  to  be 
retarded.  Although  it  is  not  very  clear  as  to  what  causes  the  retardation  at  high 
concentrations,  it  could  possibly  be  due  to  the  influence  of  interfacial  viscosities.  At  high 
concentrations  the  interface  gets  completely  saturated  and  the  intrinsic  viscosities  of  the 
surface  can  play  a  role  in  retarding  the  drop  motion.  An  effort  to  validate  this  idea 
theoretically  is  only  partially  successful. 

Finally,  Chapter  4  presents  the  analysis  of  the  displacement  flow  of  a  shear 
thinning  power-law  fluid  by  an  inviscid  fluid  in  a  Hele-Shaw  geometry.  This  problem  has 
practical  applications  in  several  methods  of  polymer  processing  such  as  injection  molding 
using  various  gases.  The  technique  of  matched  asymptotic  expansions22  is  used  to 
determine  the  flow  field  and  the  shape  of  the  interface.  The  flow  domain  is  divided  into 
three  different  regions  and  matching  conditions  between  them  provide  the  film  thickness 
left  behind  by  the  advancing  meniscus.  Unlike  the  Newtonian  case,  the  relationship 
between  capillary  number  and  film  thickness  is  not  unique  and  it  varies  from  fluid  to 
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fluid.  For  power-law  fluids,  the  depth  averaged  flow  field  is  non-Darcian.  A  linear 
stability  analysis  is  also  carried  out  to  determine  the  dispersion  relation  for  the  shear 
thinning  power-law  fluid.  Chapter  5  summarizes  the  findings  of  this  work. 


CHAPTER  2 
EXPERIMENTS 

In  this  chapter,  the  influence  of  surfactant  on  the  bubble  motion  in  a  Hele-Shaw  cell 
is  described  qualitatively,  and  the  difference  between  pressure-driven  and  buoyancy- 
driven  flows  in  the  presence  of  surfactant  is  discussed.  Subsequently,  the  results  of  the 
present  experiments  are  presented  to  support  the  arguments  described  below. 

2.1  The  Influence  of  Surfactants  on  Bubble  Motion 

As  mentioned  in  Chapter  1,  the  bubble  can  be  set  in  motion  either  by  applying  a 
pressure  gradient  in  a  horizontal  Hele-Shaw  cell  or  by  imposing  a  buoyancy  force  by 
tilting  the  cell.  When  the  bubble  is  driven  by  buoyancy  through  a  stagnant  fluid  in  a  Hele- 
Shaw  cell,  the  sweeping  motion  of  the  surrounding  liquid  in  a  bubble  fixed  coordinate 
system  is  always  from  the  front  to  the  trailing  end  of  the  bubble.  If  the  bubble  is  driven 
by  the  surrounding  liquid,  on  the  other  hand,  two  different  flow  situations  may  arise 
depending  on  the  relative  magnitude  of  the  bubble  velocity  U  to  the  average  velocity  of 
the  surrounding  fluid  V.  If  U  >  V,  the  sweeping  motion  of  the  surrounding  liquid  is 
toward  the  back  of  the  bubble  as  in  the  buoyancy  driven  flow.  Whereas  if  U  <  V,  the 
sweeping  motion  along  the  edge  of  the  bubble  is  from  the  rear  to  the  leading  end  of  the 
bubble.  In  the  absence  of  surface-active  substances,  the  surface  tension  is  constant  and 
the  two  flow  situations  may  be  indifferent.  They  are,  however,  distinctly  different  if  there 
exists  a  spatial  variation  of  surface  tension  due  to  the  presence  of  surfactant. 
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When  U  >  V,  surfactants  will  be  accumulated  at  the  trailing  end  of  the  bubble  thus 
lowering  the  surface  tension  at  the  back  of  the  bubble.  When  U  <  V,  the  flow  situation  is 
more  complicated.  At  the  edge  (or  rim)  of  the  bubble,  the  surrounding  fluid  is  convected 
from  the  back  to  the  front  of  the  bubble.  Thus,  surfactants  are  convected  toward  the 
bubble  front  along  the  edge.  In  the  thin  film  region  between  the  wall  and  the  bubble,  on 
the  other  hand,  the  dominant  convection  is  always  toward  the  back  of  the  bubble  since  the 
solid  wall  is  moving  backward  relative  to  the  bubble  and  the  liquid  film  is  very  thin. 
Therefore,  the  surfactant  distribution  is  complicated  and  may  depend  on  the  balance 
between  the  two  competing  effects.  The  studies  of  Park  and  Homsy  and  Burgess  and 
Foster  indicate  that  stagnation  points  (or  rings)  are  present  around  the  bubble  at  the 
transition  region  between  the  thin  film  region  and  the  cap  region  (i.e.,  points  s  in 
Figure  1.1).   These  stagnation  points  are  similar  to  those  in  the  front  and  the  back  of  the 
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Brethertons  bubble  in  a  capillary  tube.  Although  the  location  of  the  stagnation  points 
may  move  toward  the  cap  region  when  surfactants  are  present,  the  existence  of  the 
stagnation  points  (or  rings)  may  result  in  three  regions  of  surfactant  accumulation;  one 
near  the  front  cap  of  the  bubble  and  the  other  two  in  the  thin  film  regions  (top  and 
bottom)  near  the  back  of  the  bubble.  It  should  be  noted  that  this  complicated  flow 
situation  cannot  occur  in  a  buoyancy  driven  flow.  Therefore,  the  pressure-driven  and 
buoyancy-driven  bubble  motions  can  be  different  from  one  another  if  surfactants  are 
present  in  the  system. 

The  most  interesting  bubble  shapes  that  were  observed  by  Kopf-Sill  and  Homsy  may 
be  the  ovoids  with  the  sharper  end  pointed  in  the  flow  direction  and  the  bubbles  with  a 
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long-  or  a  short-tail.  When  the  average  liquid  velocity  V  was  low  (i.e.,  at  a  low  Ca  where 
the  surface  tension  is  dominant),  the  bubbles  were  observed  to  be  near  circular.  With 
increasing  V,  they  became  slightly  deformed  to  take  on  the  shape  of  an  ovoid.  As  V  was 
further  increased,  the  bubble  shape  went  through  a  sudden  transition  to  the  long-tail  (or 
short-tail)  shape.  Interestingly,  a  hysteresis  was  observed  in  that  the  backward  transition 
from  long-tail  to  ovoid  occurred  at  a  lower  value  of  V  than  the  forward  transition  from 
ovoid  to  long-tail  for  the  same  bubble.  One  important  aspect  which  should  be  noted  is 
that  the  ovoids  apparently  moved  with  a  velocity  smaller  than  the  average  liquid  velocity 
(i.e.,  U  <  V)  whereas  the  velocities  of  the  short-  or  long-tail  bubbles  were  larger  than  V. 
That  is,  the  sudden  increase  (or  decrease  in  case  of  a  reverse  experiment)  in  the  bubble 
velocity  cuts  across  the  U/V=l  line  in  Figure  1.3.  Thus,  the  flow  situations  for  the  ovoids 
and  the  bubbles  with  a  tail  are  different  in  that  the  direction  of  the  sweeping  motion  is 
opposite.  We  may  point  out,  however,  that  a  group  of  four  data  points  in  Figure  1 .3 
(a=0.76  cm)  does  not  follow  this  argument  as  their  relative  velocities  were  reported  to  be 
about  0.8  although  the  shapes  were  short-tail.  It  may  be  possible  that  this  discrepancy  is 
due  to  measurement  error. 

When  the  surrounding  liquid  sweeps  the  bubble  from  back  to  front  along  the  edge 
(i.e.,  when  U  <  V),  the  local  interfacial  tension  may  be  lower  in  the  front  region  of  the 
bubble  due  to  the  surfactant  accumulation  there.  Since  the  pressure  inside  the  bubble  is 
constant,  the  radius  of  curvature  at  the  bubble  front  should  be  smaller  than  that  at  the 
back  to  compensate  for  the  lower  interfacial  tension.  Therefore,  assuming  a  constant 
curvature  in  the  gap  direction  (i.e.,  the  z-direction  in  Figure  1.1),  an  ovoid  with  the 
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sharper  end  pointing  in  the  flow  direction  is  a  plausible  shape.  The  same  argument  may 
be  applicable  even  in  the  absence  of  surfactants.  In  that  case,  however,  the  difference  in 
curvature  between  the  front  and  rear  ends  of  the  bubble  will  not  be  as  large  as  the 
surfactant-influenced  case.  By  applying  the  similar  argument  to  the  case  of  U  >  V,  we 
may  also  speculate  that  a  drastic  curvature  change  in  the  rear  part  of  the  bubble  (hence 
the  long-tail  or  short-tail  bubbles)  is  possible  if  surface-active  substances  are  present. 
This  qualitative  explanation  does  not  rule  out  the  possibility  that  the  long-tail  bubbles  can 
also  exist  in  a  buoyancy  driven  flow  since  U  is  always  greater  than  V  in  that  case.  A 
probable  reason  why  it  was  not  observed  in  Maxworthy's  experiment  is  that  silicone  oil 
was  used  instead  of  glycerine-water  mixture  in  his  experiments.  As  it  was  pointed  out 
previously,  the  air-oil  interface  is  less  likely  to  be  affected  by  surface-active 
contaminants. 

Under  the  assumption  that  surface-active  contaminants  were  present,  the  overall 
behavior  of  a  bubble  observed  by  Kopf-Sill  and  Homsy  may  be  explained  as  follows. 
When  an  air  bubble  is  stagnant  in  a  contaminated  Hele-Shaw  cell,  its  entire  surface  will 
be  covered  with  surface-active  contaminants.  At  a  very  low  Ca,  the  low  shear  stress  may 
induce  only  a  small  change  in  surfactant  distribution.  Thus  the  Marangoni  effect  may  be 
global.  Although  the  Marangoni  stress  may  be  small  due  to  the  small  surface  tension 
gradient,  the  shear  stress  is  also  small  and  the  entire  bubble  surface  may  act  as  a  rigid 

24 

surface.  In  the  presence  of  the  stagnation  rings  near  the  top  and  bottom  wall  (Figure 
1.1),  three  separate  regions  of  high  surface  concentration  can  coexist  as  pointed  out 
previously.   Along  the  edge  of  the  bubble,  surfactants  may  be  swept  toward  the  bubble 


16 

front  and  consequently,  the  Marangoni  effect  may  assist  the  bubble  motion.  In  the  thin 
film  region  at  the  top  and  bottom  of  the  walls,  on  the  other  hand,  surfactants  are  swept 
backward  retarding  the  bubble  motion.  When  Ca  is  small,  the  viscous  stress  along  the 
bubble  edge  may  not  be  as  significant  as  that  in  the  thin  film  region  which  is  inversely 
proportional  to  the  film  thickness.  Therefore  the  drag  in  the  wetting  film  is  dominant 
retarding  the  bubble  motion  significantly.  Consequently,  the  relative  bubble  velocity 
(U/V)  may  remain  much  smaller  than  1 .  Although  the  bubble  motion  can  be  severely 
retarded,  the  shape  change  may  be  insignificant  due  to  the  small  spatial  variation  of  the 
surface  tension.  Consequently,  the  bubble  shape  may  be  near-circular. 

As  Ca  is  increased,  the  larger  viscous  stress  along  the  bubble  edge  increases  the 
surface  tension  gradient  resulting  in  larger  Marangoni  stress.  In  the  thin  film  region,  on 
the  other  hand,  the  relative  increase  of  the  drag  may  not  be  as  significant  since  the  film 
thickness  increases  with  the  bubble  velocity.  Therefore,  the  relative  bubble  velocity  may 
increase  with  increasing  Ca  although  its  rate  of  increase  may  not  be  very  large. 
Furthermore,  the  surfactant  concentration  at  the  leading  end  of  the  bubble  will  also 
increase  with  Ca.  Consequently,  the  surface  tension  will  become  lower  at  the  leading  end 
creating  a  larger  curvature.  Therefore,  the  ovoid  with  the  sharper  end  pointing  in  the 
direction  of  the  flow  may  be  a  favorable  shape. 

With  increasing  Ca,  the  film  thickness  increases  and  the  location  of  the  stagnation 
rings  may  move  toward  the  tip  of  the  bubble  edge.  Eventually  the  two  stagnation  rings 
may  collapse  into  one  along  the  tip  of  the  bubble  edge.  Then  the  three  spearate  regions 
of  high  surfactant  concentration  cannot  coexist.  Consequently,  the  large  drag  in  the  thin 
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film  region  may  decrease  significantly  inducing  higher  bubble  velocity  which  may  be 
greater  than  the  average  velocity  of  the  surrouding  fluid  (i.e.,  U/V  >  1).  Once  the  relative 
bubble  velocity  becomes  greater  than  1 ,  all  surfactants  will  be  collected  at  the  trailing  end 
of  the  bubble  lowering  the  surface  tension  there.  Consequently,  a  larger  curvature  is 
required  forming  a  high  curvature  tail.  The  size  (or  length)  of  the  tail  may  vary 
depending  upon  the  surfactant  distribution  at  the  rear  end.  In  an  extreme  case  where  the 
bubble  is  large  enough  and  the  surfactant  transport  is  limited  by  bulk  diffusion,  the 
leading  end  of  the  bubble  may  be  virtually  free  from  surfactant  due  to  convection  whereas 

the  trailing  end  is  full  of  surfactants.  '  In  such  a  case,  a  long  cylindrical  tail  may  be 
conceivable  as  a  form  of  minimizing  viscous  dissipation  while  maintaining  a  large 
curvature  for  the  pressure  balance  across  the  interface.  The  Marangoni  effect  is  greater  if 
the  surface  curvature  is  larger.  Consequently,  its  strong  influence  may  persist  to  a  higher 
value  of  Ca  if  the  bubble  is  smaller.  Therefore,  the  critical  value  of  Ca  where  the 
transition  occurs  may  increase  with  decreasing  bubble  size. 

If  the  capillary  number  is  decreased  gradually  from  a  large  value  to  a  small  value,  the 
reverse  trend  may  be  observed  in  which  the  bubble  goes  through  a  transition  from  a  long- 
or  short-tail  bubble  to  an  ovoid  and  eventually  to  a  near-circle.  The  forward  and 
backward  transitions,  however,  are  not  necessarily  reversible  as  the  adsorption  and 
desorption  kinetics  of  surfactants  may  not  be  reversible.  Thus,  the  hysteresis  in  the 
transitional  behavior  can  occur  if  such  a  change  is  driven  by  the  presence  of  surfactants. 
The  argument  presented  above  is  rather  simple  and  only  qualitative.  Nevertheless,  it 
apparently  accounts  for  most  of  the  observations  of  Kopf-Sill  and  Homsy.     In  the 
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following  section,  the  results  of  the  present  experiments  are  given  as  a  supporting 
evidence. 

2.2  Experiments 

The  Hele-Shaw  cell  consisted  of  two  1/2-in  thick  Pyrex  glass  plates  separated  by  a 
rubber  gasket  of  0.9  mm  or  1.8  mm  in  thickness.  The  effective  cell  dimension  was  17.8 
cm  by  86.4  cm.  A  silicone  oil  with  the  measured  viscosity  and  surface  tension  of  97  cp 
and  21  dyn/cm  was  used  as  the  driving  fluid.  The  cell  had  a  deep  groove  channel  at  both 
ends  so  that  a  uniform  distribution  of  the  oil  along  the  transverse  direction  could  be 
ensured  when  it  was  injected.  The  cell  also  had  an  injection  port  at  one  end  of  the  top 
plate  so  that  an  air  bubble  or  a  water  drop  could  be  introduced  into  the  cell  using  a 
syringe  without  creating  extra  bubbles  or  drops.  The  planform  diameter  of  the  stagnant 
bubbles  and  drops  was  controlled  to  be  at  about  1.3  cm  and  2.1  cm  throughout  the 
experimental  study.  These  sizes  were  chosen  to  be  about  the  same  as  those  of  Kopf-Sill 
and  Homsy  for  direct  comparison  of  the  results.  While  the  size  variation  among  the  small 
bubbles  was  small  as  they  ranged  between  1 .3  cm  and  1 .4  cm,  the  size  control  was  more 
difficult  with  larger  bubbles  and  their  size  ranged  between  2.2  cm  to  2.5  cm. 

Water  drops  containing  an  organic  surfactant  at  various  levels  of  concentration  were 
prepared  by  dissolving  sodium  dodecyl  sulfate  (SDS)  in  distilled  water.  Since  the 
surfactant  influence  is  more  significant  at  low  concentrations,  the  SDS  concentration  was 
set  at  5%,  10%,  20%,  and  100%  of  the  critical  micelle  concentration  (CMC),  respectively. 
The  CMC  of  SDS  in  water  at  25°C  is  8.2  mmol/l  which  is  equivalent  to  0.236%  by 
weight.    The  interfacial  tension  measured  at  room  temperature  by  the  Wilhelmy  plate 
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method  is  listed  in  Table  2.1  at  various  SDS  concentrations.  As  the  present  data  for  air- 
water  interface  is  found  to  be  consistent  with  those  reported  in  literature,  the  data  for  the 
silicone  oil-water  interface  may  also  be  reliable.  It  can  be  deduced  from  the  data  that  the 
rate  of  change  of  interfacial  tension  is  greater  at  a  lower  concentration  indicating  a 
stronger  Marangoni  effect  at  a  lower  concentration. 

Once  an  air  bubble  (or  a  water  drop)  of  predetermined  size  was  positioned  at  one  end 
of  the  cell,  the  silicone  oil  was  driven  by  a  peristaltic  pump  at  various  flow  rates  to  induce 
the  bubble  motion.  The  translational  velocity  and  the  shape  of  the  bubble  (or  drop)  were 
then  recorded  for  each  value  of  oil  flow  rate.  For  the  prescribed  dimensions,  the  values  of 
both  a/L  and  h/a  were  smaller  than  0.14.  Thus  the  influence  of  the  side  wall  (i.e.,  the 
rubber  gasket)  may  be  small  and  the  depth-averaged  Hele-Shaw  equations  may  be 
applicable  to  describe  the  bubble  motion  if  the  system  is  free  from  surfactants.  In  order 
to  avoid  an  excessive  pressure  drop  at  high  flow  rates,  most  of  the  experiments  were 
conducted  with  the  thick  gasket  (1.8  mm)  and  only  a  part  of  the  air  bubble  experiments 
used  the  thin  gasket  (0.9  mm)  for  comparison. 

The  Bond  number  (Bo)  which  is  defined  as  pgh^/a  is  a  measure  of  the  relative 
importance  of  gravity  to  surface  tension,  and  a  large  Bo  may  indicate  asymmetry  of  the 
bubble  in  the  gap.  In  the  experiments  of  Kopf-Sill  and  Homsy,  Bo  was  as  small  as  0. 1 
whereas  it  was  about  1 .5  in  Maxworthy's.  In  the  present  experiment,  it  was  either  0.4  for 
the  0.9  mm  thick  gasket  or  1.5  for  the  1.8  mm  thick  gasket.  The  large  Bo  was  caused  not 
only  by  the  thick  gasket  but  also  by  the  low  interfacial  tension.  While  the  gravity  effect 
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might  not  be  significant  when  Bo  =  0.4,  uncertainty  existed  when  Bo  was  as  large  as  1.5. 
The  experimental  results  for  air  bubbles,  however,  were  indifferent  in  terms  of  both 
bubble  velocity  and  shape  to  whether  Bo  was  0.4  or  1.5.  Thus,  the  gravity  effect  might 
not  have  significant  influence  even  when  Bo  =  1.5.  In  reporting  the  results  in  the 
following  section,  no  distinction  has  been  made  between  the  data  for  Bo  =  1 .5  and  0.4  in 
order  to  avoid  unnecessary  complication. 

The  size  of  air  bubbles  changed  slightly  during  the  experiment  due  to  the  pressure 
drop  in  the  cell.  Kopf-Sill  and  Homsy  observed  similar  behavior  and  reported  the  volume 
change  of  air  bubbles  up  to  20%.  In  the  present  experiment,  not  only  was  the  gap 
thickness  of  the  cell  slightly  larger  but  also  the  liquid  viscosity  was  lower.  Furthermore, 

the  liquid  flow  rate  was  0.75  cm^/sec  when  the  capillary  number  was  as  high  as  110. 
Consequently,  the  pressure  drop  was  smaller  and  the  change  in  the  bubble  volume  was 
less  than  about  5%.  Despite  the  change  in  the  bubble  volume,  variation  in  the  bubble 
velocity  was  not  detectable.  In  case  of  water  drops,  neither  the  volume  nor  the  velocity 
changed. 

The  description  in  Section  2.1  regarding  the  surfactant  influence  on  the  bubble 
motion  assumes  the  existence  of  a  thin  oil  film  between  the  wall  and  the  moving  bubble 
or  drop.  When  oil  is  the  wetting  phase,  the  film  thickness  is  uniform  in  the  direction  of 
bubble  or  drop  motion  (i.e.,  x  direction  in  Figure  1.1)  whereas  it  varies  in  the  cross-stream 
direction  (i.e.,  v  direction  in  Figure  1.1).  Therefore,  it  may  be  possible  in  principle  to 
prove  the  existence  of  the  thin  film  by  an  optical  method  which  creates  interference 
fringes.  Our  attempt  to  observe  the  interference  fringes,  however,  was  unsuccessful  due 
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to  various  difficulties  involved  with  the  large  size  of  the  cell.     Thus  we  provide  only  a 
qualitative  evidence  to  support  the  assumption. 

A  photograph  converted  from  a  video  image  is  shown  in  Figure  2.1  for  a  water  drop 
of  about  2  mm  in  diameter  which  is  sitting  on  a  Pyrex  glass  plate  immersed  in  the 
silicone  oil.  In  order  to  improve  contrast,  a  blue  food  coloring  was  added  to  the  water 
drop.  Due  to  the  difficulty  in  obtaining  perfectly  horizontal  alignment  for  the  contact 
angle  measurement,  the  water  drop  was  viewed  with  a  small  positive  angle  from  the 
horizontal  and  a  shadow  of  the  drop  was  created  in  front  by  indirect  back  lighting.  The 
lower  part  of  the  picture  is  the  shadow  of  the  drop  and  the  bright  line  between  the  drop 
and  its  shadow  represents  the  three-phase  contact  line.  At  the  edge  of  the  contact  line,  the 
apparent  static  contact  angle  is  shown  to  be  greater  than  90°  when  measured  from  the 
water-drop  side.  This  may  suggest  that  silicone  oil  is  the  wetting  phase  under  dynamic 
conditions. 

When  a  water  drop  is  stationary  in  the  Hele-Shaw  cell  at  the  beginning  of  each 
experiment,  the  water  drop  may  have  direct  contact  with  the  Pyrex  glass  wall  creating  a 
three-phase  contact  line.  As  we  induced  the  drop  motion  by  injecting  the  silicone  oil  into 
the  cell,  the  trailing  edge  of  the  drop  exhibited  briefly  a  ragged  contact  line  pinned  on  the 
wall.  The  contact  line  soon  detached  from  the  wall  and  the  water  drop  took  on  a  steady 
shape  with  a  perfectly  smooth  edge  around  the  entire  drop  (Figure  2.2).  If  a  moving 
contact  line  existed,  its  motion  might  have  been  jerky  at  times  unless  the  solid  surface 
was  perfectly  homogeneous.  In  the  present  experiments,  the  drop  motion  was  always 
smooth  excluding  the  brief  moment  at  the  start  of  motion.  In  case  of  air  bubbles,  even  the 
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Figure  2.1     A  water  drop  sitting  on  a  Pyrex  glass  plate  immersed  in  silicone  oil. 
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start-up  motion  was  smooth  without  any  evidence  of  contact  line.  These  observations 
along  with  the  observation  on  the  apparent  static  contact  angle  seem  to  support  the 
assumption  regarding  the  existence  of  thin  oil  films. 

2.3  Results  and  Discussion 

In  Figure  2.2,  shapes  of  a  water  drop  containing  surfactant  are  given  sequentially 
with  increasing  capillary  number.  These  pictures  were  converted  from  the  video  images 
of  a  water  drop  which  was  2.5  cm  in  diameter  and  contained  sodium  dodecyl  sulfate 
(SDS)  at  20%  of  the  critical  micelle  concentration  (CMC).  The  direction  of  the  flow  is 
downward  from  top  to  bottom  of  the  figure.  In  order  to  improve  contrast  a  blue  food 
coloring  was  added  to  the  drop.  It  was  confirmed  by  measurements  that  the  food  coloring 
did  not  affect  the  surface  tension.  All  shapes  were  observed  to  be  steady  excluding  (2.2e) 
which  was  a  transitional  shape  from  (2.2d)  to  (2.2f).  Figure  2.2a  is  the  image  of  a 
stationary  drop  at  Ca  =  0.  Due  to  the  typical  distortion  of  a  video  image  in  the  vertical 
direction,  this  circular  drop  appears  to  be  slightly  elliptic.  As  Ca  was  increased,  the  drop 
was  elongated  sideways  (2.2b).  Further  increase  in  Ca  resulted  in  further  stretching  of 
the  drop  in  the  transverse  direction  (2.2c).  This  shape  is  apparently  asymmetric  about  the 
vertical  axis  (or  the  x-axis  in  Figure  1.1).  Nevertheless,  its  dimensions  remained 
unchanged  and  moved  with  a  steady  velocity.  Hence  a  steady  shape. 

It  should  be  pointed  out  that  the  severely  stretched  drop  shape  (2.2c)  occurred 
only  when  the  drop  was  larger  than  about  2.1  cm  in  diameter  and  when  it  contained  SDS 
at  a  concentration  between  5%  and  20%  of  CMC.  It  did  not  occur  with  the  pure  water 
drop  nor  with  the  one  containing  SDS  at  100%  CMC  although  their  sizes  were  larger  than 
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(a) 


(b) 


(c) 


(d) 


Figure  2.2  Evolution  of  water  drop  shape  with  increasing  capillary  number,  (a)  Near 
circular  drop;  (b)  flattened  drop;  (c)  severely  stretched  drop;  (d)  ovoid;  (e) 
transitional  shape;  and  (f)  short  tail  drop 
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2.1  centimeters.  Furthermore,  when  Ca  was  increased  further  from  the  value  for  Figure 
2.2c,  the  drop  broke  into  two  instead  of  showing  a  transition  to  an  ovoid  or  short-tail.  In 
order  to  obtain  the  steady  shapes  at  a  higher  Ca  (i.e.,  Figure  2.2d  or  2.2f)  with  the  large 
drops  containing  SDS  at  5%  to  20%  CMC,  the  Ca  had  to  be  set  to  a  large  value  from  the 
beginning.  Then  the  drop  evolved  directly  from  (2.2a)  to  (2.2d)  or  (2.2f)  bypassing  the 
severe  sideways  stretching  (2.2c).  In  case  of  small  drops  (1.3  or  1.35  cm  in  diameter)  or 
the  large  drops  with  0%  or  100%  CMC,  on  the  other  hand,  the  severely  stretched  shape 
(2.2c)  did  not  occur  and  the  evolution  of  steady  shapes  was  smooth  from  (2.2a)  to  (2.2f) 
with  the  gradual  increase  of  Ca.  It  may  be  interesting  to  note  that  the  drop  shape  (2.2c) 
has  a  negative  curvature  at  the  leading  edge  as  in  the  case  of  "Tanveer"  bubbles  although 
the  trailing  edge  also  shows  a  negative  curvature. 

The  steady  shape  (2.2d)  which  occurred  at  a  higher  Ca  than  (2.2b)  or  (2.2c)  is 
different  from  others  as  the  curvature  of  its  leading  edge  (i.e.,  the  lower  part  of  the  drop  in 
Figure  2.2d)  is  apparently  larger  than  that  of  the  trailing  edge.  We  may  call  this  shape 
"ovoid"  with  its  sharper  end  pointing  in  the  direction  of  flow  which  is  equivalent  to  the 
sketch  given  in  Figure  1.2d.  The  relative  velocity  (i.e.,  U/V)  of  this  drop  was  still  smaller 
than  1 .  When  Ca  was  further  increased,  the  ovoid  went  through  a  rapid  transition  to  a 
short-tail  (2.2f)  which  moved  faster  than  the  average  velocity  of  the  surrounding  fluid 
(i.e.,  U/V  >  1).  This  transition  always  occurred  through  a  transitional  shape  given  in 
Figure  2.2e. 

In  Figure  2.3,  2.4,  and  2.5,  the  bubble  or  drop  velocity  relative  to  the  average 
velocity  of  the  surrounding  fluid  is  given  as  a  function  of  capillary  number  Ca.    Also 
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Figure  2.3  U/V  vs  Capillary  number  for  air  bubbles. 
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specified  in  the  figures  are  the  steady  shapes  of  the  bubbles  or  drops  classified  according 
to  Figure  2.2.  The  capillary  number  in  these  plots  is  defined  as  uV/a  where  fa  and  V  are 
the  viscosity  and  the  average  velocity  of  the  silicone  oil  respectively  and  a  the 
equilibrium  interfacial  tension.  In  case  of  air  bubbles  (Figure  2.3),  a  is  21  dyn/cm 
whereas  a  for  the  water  drops  (Figures  2.4  and  2.5)  varies  depending  on  the  SDS 
concentration.  The  equilibrium  values  at  each  SDS  concentrations  are  given  in  Table  2.1 
and  these  respective  values  have  been  used  in  calculating  the  capillary  number. 

The  relative  velocities  of  air  bubbles  described  in  Figure  2.3  are  for  two  different 
sizes  (a=1.27  cm  and  2.54  cm).  Overlaid  in  the  figure  are  the  data  taken  from  Kopf-Sill 
and  Homsy  for  similar  size  bubbles.  In  their  experiments,  both  small  and  large  bubbles 
experienced  a  sudden  increase  in  the  relative  velocity  accompanied  by  the  shape  change. 
The  transition  for  the  small  bubble  occurred  when  Ca  was  about  10"  whereas  it  occurred 
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at  Ca  w  5  X  10  for  the  large  bubble.  Such  transitions,  however,  were  not  observed  in  the 
present  experiment  which  covered  the  range  of  Ca  between  2  X  10"3  and  1.1  X  10"2.  The 
relative  velocity  was  always  larger  than  one  for  both  bubbles,  although  the  larger  bubble 
was  moving  slightly  faster  than  the  smaller  one.  Furthermore,  the  bubble  shapes  were 
always  near  circular  unlike  the  experiments  of  Kopf-Sill  and  Homsy. 

The  lower  limit  of  Ca  (2  X  10" )  was  set  by  the  sensitivity  of  the  pump  whereas  the 
upper  limit  (1.1  X  10  )  was  determined  by  the  maximum  flow  rate  which  was  imposed 
by  the  pump  capacity  or  the  upper  limit  of  pressure  drop  in  the  cell.  The  lower  limit  is 
well  below  the  critical  Ca  of  Kopf-Sill  and  Homsy  at  which  the  small  bubble  experienced 
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the  velocity/shape  transition.  The  critical  value  for  the  large  bubble  (Ca  «  5  X  10  ),  on 
the  other  hand,  was  smaller  than  the  lower  limit.  Thus,  the  present  experiment  could  not 
have  detected  the  transition  if  it  had  occurred  at  a  very  low  Ca.  It  should  be  noted, 
however,  that  in  the  experiments  of  Kopf-Sill  and  Homsy,  the  velocity  transition  was 
always  accompanied  by  the  shape  change  and  the  bubbles  with  the  relative  velocity  (U/V) 
greater  than  one  had  either  a  long-  or  a  short-tail.  In  the  present  experiments,  on  the 
other  hand,  the  bubble  shapes  were  near-circular  although  their  relative  velocities  (U/V) 
were  greater  than  about  1.5.  Thus,  it  is  unlikely  that  these  near-circular  bubbles  had 
experienced  the  sudden  transition  at  a  smaller  Ca  than  the  lower  limit  or  they  would 
experience  it  at  a  higher  Ca  than  the  upper  limit. 

Considering  the  slight  deviation  from  a  circular  shape  and  the  viscous  dissipation  at 

the  bubble  edges  which  was  pointed  out  by  Maxworthy,  the  bubble  velocities  in  the 
present  experiments  appear  to  be  in  reasonable  agreement  with  the  prediction  of  Taylor- 
Saffman.  The  distinct  differences  between  the  results  of  the  present  experiments  and  of 
Kopf-Sill  and  Homsy  are  apparently  in  accordance  with  the  argument  described  in  section 
2.1  that  they  were  probably  caused  by  the  influence  of  surface-active  contaminants. 

In  Figures  2.4  and  2.5,  the  results  for  water  drops  are  given  for  two  different  sizes. 
In  these  experiments  the  1.8  mm  gasket  was  used  to  minimize  the  pressure  in  the  cell  at 
high  flow  rate,  and  the  maximum  attainable  capillary  number  was  determined  by  the 
capacity  of  the  pump  (Cole-Parmer,  Masterflex  Model  7520-25).  Due  to  the  variation  of 
the  equilibrium  interfacial  tension  with  the  SDS  concentration,  the  maximum  capillary 
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Table  2.1    Variation   of  the   interfacial   tension   with   sodium   dodecyl    sulphate 
concentration  at  25°C 


SDS  concentration  in  water 


(%ofCMC) 


Interfacial  Tension  (dyn/cm) 
air-water  Silicone  oil-water 


0 

5 

10 
20 
100 


72.8 
58.7 
53.8 
41.7 
33.0 


38.8 

34.3 

32.1 

25.9 

13.7 

30 


number  also  varied  accordingly  for  the  given  maximum  flow  rate.  In  case  of  the  100% 
CMC  water  drop,  the  equilibrium  interfacial  tension  was  as  low  as  13.7  dyn/cm  whereas 
that  of  pure  water  and  oil  was  38.8  dyn/cm  (Table  2.1).  Consequently,  the  maximum 

capillary  number  for  the  former  (3.7  X  10  )  was  higher  than  that  for  the  latter  (1 .3  X  10  ) 
by  a  factor  of  about  2.8. 

The  relative  velocities  of  the  small  water  drops  (a  =  1.30  cm  or  1.35  cm)  are  shown 
in  Figure  2.4.  The  labels  (a)-(f)  indicate  the  drop  shapes  as  indicated  in  Figure  2.2. 
Unlike  the  air  bubbles,  they  were  moving  very  slowly  at  low  capillary  numbers  as  in  the 
observations  of  Kopf-Sill  and  Homsy.  It  is  interesting  to  note  that  the  water  drops  which 
were  supposed  to  be  pure  were  also  moving  very  slowly  at  low  Ca  although  they  were 
expected  to  be  as  fast  as  the  air  bubbles  in  Figure  2.3.  A  plausible  explanation  for  this 
behavior  may  be  that  the  water  drops  still  contained  some  surfactant  despite  the 
precaution  we  had  taken  to  prevent  contamination.  When  the  SDS  concentration  was 
between  5%  and  20%  of  CMC,  the  relative  velocities  were  smaller  than  one  even  at  the 
upper  limit  of  Ca.  Furthermore,  the  drop  velocities  were  indistinguishable  and  the  shapes 
were  near-circular  or  elliptic  (Figures  2.2a  and  2.2b).  Although  it  was  expected  that  the 
relative  velocities  of  these  drops  would  increase  with  Ca  and  become  eventually  larger 
than  one,  the  transition  could  not  be  observed  as  it  might  occur  at  a  higher  Ca  than  the 
attainable  upper  limit. 

In  case  of  the  "pure"  drop  or  the  drop  with  SDS  at  100%  CMC,  on  the  other  hand, 
the  transition  could  be  observed  in  that  the  drops  took  on  the  short-tail  shape  (Figure  2.2f) 
and  moved  with  a  relative  velocity  greater  than  one.  Unlike  the  observations  of  Kopf-Sill 
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Figure  2.  4         UA''  vs  Capillary  number  for  small  water  drops. 
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and  Homsy,  however,  the  change  (or  the  transition)  was  gradual  rather  than  abrupt. 
Although  it  is  not  clear  what  caused  this  difference  in  the  transitional  behavior,  it  may  be 
partly  due  to  the  differences  between  the  two  experiments.  In  the  current  experiment  of 
water  drop  and  silicone  oil,  the  surfactant  was  present  inside  the  water  drop  whereas  in 
case  of  Kopf-Sill  and  Homsy's  experiment  of  air  bubble  and  aqueous  solution,  the 
surfactants  (or  surface-active  contaminants)  might  be  present  in  the  aqueous  solution 
outside  the  bubble.  Therefore,  the  mechanism  for  surfactant  transport  to  the  bubble  or 
drop  surface  could  be  different.  Furthermore,  considering  the  differences  in  the  fluid 
properties,  quantitative  match  between  the  two  experiments  may  not  be  expected.  These 
differences  could  be  also  the  reason  why  long-tail  drops  were  not  observed  in  the  present 
study. 

The  large  water  drops  showed  the  similar  trend  as  the  small  ones  in  that  the  relative 
velocities  were  smaller  than  one  at  low  Ca  (Figure  2.5).  Unlike  the  case  of  small  drops, 
however,  the  transitional  behavior  is  evident  for  all  water  drops  in  Figure  2.5.  This 
observation  is  consistent  with  the  trend  in  Figure  1.3  in  that  the  critical  Ca  for  the 
velocity/shape  transition  decreased  with  increasing  bubble  size.  As  we  mentioned  in 
section  2.1,  the  surfactant  influence  (or  the  Marangoni  effect)  should  be  greater  if  the 
surface  curvature  is  larger.  Consequently,  the  strong  influence  of  surfactant  may  persist 
to  a  higher  value  of  Ca  if  the  drop  is  smaller.  Apparently,  the  critical  Ca  for  the  large 
drops  was  smaller  than  the  upper  limit  of  Ca  of  the  present  experiments  whereas  that  for 
the  small  ones  was  larger  than  the  upper  limit.  The  variation  in  the  drop  shape  was  also 
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consistent  in  that  it  was  near-circular,  elliptic  or  ovoid  when  U/V  <  1,  whereas  all  drops 
became  short-tailed  when  U/V  >  1.  Figure  2.5  also  indicates  that  the  critical  Ca  for  the 
velocity  and  shape  transition  increases  with  increasing  SDS  concentration  with  the 
exception  of  the  drop  at  100%  CMC.  Although  only  a  quantitative  analysis  will  provide 
proper  answer  to  such  a  trend,  a  qualitative  explanation  may  be  given  as  follows.  For 
most  organic  surfactants  the  adsorption  kinetics  is  fast.  Therefore,  if  the  surfactant 
concentration  is  low,  its  transport  may  be  limited  by  bulk  diffusion.  At  a  reasonably  high 
Ca,  the  surfactant  distribution  may  be  localized  due  to  convection  and  consequently,  the 
Marangoni  effect  is  also  restricted  to  the  local  region.  As  long  as  the  surface 
concentration  gradient  is  smaller  than  the  maximum  allowable  value  before  the 
monolayer  collapse  for  the  given  surfactant,  the  area  covered  with  the  surfactant  may 
grow  with  increasing  bulk  concentration.  Consequently,  the  Marangoni  effect  may  also 
increase  with  the  bulk  concentration.  If  the  bulk  concentration  is  higher  than  a  certain 
limit,  however,  surfactant  transport  may  not  be  limited  by  bulk  diffusion  any  more  and 
the  surface  concentration  gradient  will  start  to  decrease.  Thus  the  surfactant  influence 
will  also  diminish  with  further  increase  in  bulk  concentration.  It  seems  that  the  critical 
bulk  concentration  for  the  present  experiment  at  which  the  surfactant  influence  starts  to 
decrease  may  be  between  20%  and  100%  of  CMC.  In  an  extreme  case,  if  the  bulk 
concentration  is  very  high,  the  surface  may  be  remobilized  due  to  the  lack  of  surface 
concentration  gradient.  We  have  also  conducted  an  experiment  with  a  water  drop  which 
contained  SDS  at  10  times  the  CMC  to  investigate  whether  the  remobilization  would 
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Figure  2.6  U/V  vs  Capillary  number  for  Fluorocarbon  oil  drops. 
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occur  at  that  concentration.    The  result,  however,  almost  coincided  with  that  for  1 00% 
CMC  drop  in  Figure  2.5  in  that  the  relative  drop  velocity  was  still  very  low  at  small  Ca. 

Finally,  an  experiment  is  conducted  using  an  oil-oil  system.  Although,  surface  active 
substances  are  present  as  contaminants  in  this  system,  they  are  not  expected  to  adsorb  at 
the  interface  hence  their  influence  on  the  drop  motion  will  be  insignificant.  Fluorocarbon 
oil  is  chosen  because  it  is  immiscible  with  silicone  oil.  The  interfacial  tension  was 
measured  by  Wilhelmy  plate  method  and  is  found  to  be  close  to  5  dyns/cm.  Pressure 
gradient  is  imposed  to  study  the  motion  of  a  Fluorocarbon  oil  drop  placed  in  Hele-Shaw 
cell  filled  with  silicone  oil.  Figure  2.6  presents  the  normalised  velocity  of  fluorocarbon 
oil  drop  as  a  function  of  capillary  number.  A  nominal  value  of  5  dyns/cm  is  used  in 
calculating  the  capillary  number.  These  results  indicate  that  the  normalised  velocity  is 
always  greater  than  one  and  the  drop  velocity  is  close  to  the  predictions  of  Taylor- 
Saffman  theory.  Unlike  the  water  drops,  the  shape  of  the  fluorocarbon  oil  drop  is  always 
found  to  be  circular.  These  results  reinforce  the  idea  that  the  retardation  experienced  by 
drops  or  bubbles  in  an  aqueous  environment  is  due  to  the  influence  of  surfactants. 

2.4  Summary  and  Conclusions 

Under  the  assumption  that  surface-active  contaminants  were  the  primary  reason  for 
the  interesting  behaviors  of  the  bubbles  in  a  Hele-Shaw  cell  observed  by  Kopf-Sill  and 
Homsy,  similar  experiments  were  conducted  using  the  fluid  combinations  of  air-silicone 
oil  and  water-silicone  oil.  These  systems  were  chosen  to  delineate  the  surfactant 
influence  systematically.  In  case  of  air  bubbles  in  silicone  oil  for  which  the  surfactant 
influence  was  likely  to  be  insignificant,  the  relative  velocities  were  observed  to  be  as  high 
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as  the  Taylor-Saffman  prediction.  Furthermore,  the  unusual  bubble  shapes  such  as  long- 
or  short-tail  bubbles  were  not  observed.  In  case  of  water  drops  which  contained 
predetermined  amount  of  an  organic  surfactant  (sodium  dodecyl  sulfate),  on  the  other 
hand,  very  low  translational  velocities  as  well  as  the  unusual  shapes  of  Kopf-Sill  and 
Homsy  were  observed. 

Although  it  is  impossible  to  confirm  the  presence  of  surfactants  in  the  past 
experiments,  the  present  results  are  consistent  with  those  of  Kopf-sill  and  Homsy  and  it 
appears  that  most  of  the  observations  by  them  may  be  due  to  the  influence  of  surface- 
active  contaminants. 


CHAPTER  3 
ESTIMATION  OF  BUBBLE  VELOCITY 

3.1  Introduction 

This  chapter  presents  a  theoretical  calculation  for  the  translational  velocity  of  a 
bubble  which  is  retarded  by  the  surfactant  influence.  In  order  to  simplify  the  analysis,  the 
bubble  shape  is  assumed  a  priori  to  be  elliptic.  It  is  also  assumed  that  the  surrounding 
fluid  wets  the  solid  surface  thereby  forming  a  thin  film  between  the  bubble  and  the  plates 
of  the  Hele-Shaw  cell.  For  a  relatively  small  bubble  (i.e.,  b/L  «  1)  yet  large  enough  to 
neglect  the  edge  effect  (h/b  «  1),  the  entire  bubble  surface  tends  to  be  rigid  due  to  the 
surface  active  substance.  Consequently,  a  large  drag  arises  in  the  thin  film  region  thus 
retarding  the  bubble  motion  significantly.  The  calculated  velocity  of  the  bubble  is 
apparently  in  reasonable  agreement  with  our  own  experimental  observations  and  also 
with  those  of  Kopf-Sill  and  Homsy.6  In  addition,  the  present  calculation  provides  an 
explanation  for  the  evolution  of  the  bubble  shape  with  increasing  capillary  number. 

3.2  Theoretical  Analysis 

A  bubble  or  drop  in  a  Hele-Shaw  cell  assumes  a  circular  plan  form  when  it  is  stagnant 
because  of  the  interfacial  tension.  Once  the  bubble  is  set  in  motion,  its  shape  deviates 
from  the  circular  plan  form  and  takes  on  a  steady  shape  depending  on  the  flow  conditions 
(Figure  1.1).  In  the  absence  of  surface  tension  effect,  the  analysis  by  Taylor  and 
Saffman1  predicts  an  elliptic  shape  if  the  bubble  size  is  much  smaller  than  the  width  of 
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the  cell  (i.e.,  b/L  «  1).  If  surface  active  substances  are  present,  the  bubbles  are  not 
necessarily  elliptical  and  may  take  on  various  interesting  shapes  depending  on  the  flow 
conditions6.  Nevertheless,  we  assume  an  elliptic  plan  form  since  the  shape  distortion  is 
small  at  low  capillary  numbers.  While  the  bubble  shape  should  be  determined  as  a  part  of 
solution  in  a  rigorous  analysis,  the  present  analysis  is  an  approximation  in  which  the 
bubble  shape  is  assumed  a  priori.  Unlike  a  full  three-dimensional  analysis  that  may  be 
possible  only  numerically,  the  present  analysis  provides  an  analytic  description  for  the 
translational  velocity  of  the  bubble. 

As  indicated  in  Figure  1 . 1 ,  the  two  principal  axes  of  the  elliptic  bubble  are  assumed  to 
be  aligned  with  the  flow  and  the  transverse  directions.  The  surrounding  fluid  wets  the 
solid  surface  thus  forming  a  thin  liquid  film  between  the  plates  and  the  bubble.  In  the 
presence  of  surface  active  substances,  the  Marangoni  effect  resulting  from  the  surface 
tension  gradient  may  complicate  the  flow  field  near  the  bubble.  When  the  bubble  is 
small,  however,  it  may  be  simplified  since  the  entire  bubble  surface  may  become 
rigid19,20.  Here  we  consider  the  case  of  a  small  bubble  in  which  the  entire  bubble  surface 
is  assumed  to  be  rigid. 

The  flow  field  slightly  away  from  the  bubble  is  known  to  be  parabolic  in  the  xz-plane 
(i.e.,  across  the  gap).    Thus,  the  velocity  field  can  be  given  as1,7 


v  =  1.5v 

avg 


V       h2 ) 


Here  vavg  is  the  gap  (or  depth)-averaged  velocity  field  that  satisfies  the  following  Hele- 
Shaw  equations: 
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V-vavg  =  0  3.2 


12u 
h2  ' 


Vp  =  --fvavg  3.3 


Here  p  is  the  gap-averaged  pressure  field,  p.  is  the  fluid  viscosity  and  h  is  the  gap  of  the 
Hele-Shaw  cell.  It  is  apparent  from  Eqns.  (3.2)  and  (3.3),  that  the  average  flow  field  in  the 
xy-plane  can  be  described  as  a  potential  flow.  In  an  elliptic  cylindrical  coordinate  system 
depicted  in  Figure  3.1,  the  complex  potential  Q  for  the  flow  past  an  elliptic  cylinder  in  a 
bubble  fixed  frame  of  reference  is  given  by16 

Q  =  <D  +  ,"F=    -b(y-u)£±^J?—x\j?!__JLL]  3.4 

v  '      2  [jfc  +  1      Jt-lJ 

where  y  =  %  +  i  r|  and  k  the  shape  factor  defined  as  k  =  a/b.  V  and  U  are  the  average 

velocity  of  the  surrounding  fluid  and  the  bubble,  respectively.  From  the  gradient  of  the 

velocity  potential  O  or  the  stream  function  x¥,  the  velocity  components  in  \  and  r| 

directions  are  given  as 


k+ 1 1   e^        e~^  1  Uc 

v^g  =  c(V-U) 1 costi  +  — sinh^cosri 

^  V  '  2q  [k  +  \     k-\]         '      q  S        ' 

k  +  \[  e^        e~^  1  Uc 

wavg  =  -c(V-  U) 1 h isinri-  —  coshtsinri 

1  V  '  2q  [k  +  \     k-\\        '      q  h        ' 


3.5 


3.6 


Here  c  =  J|a    -  b  |  and  q  is  the  scale  factor  for  the  elliptic  cylindrical  coordinate  system 

defined  as  q2  -  c1  (sinh2^+sin2r|).  As  it  was  pointed  out  previously,  Eqns  (3.5)  and  (3.6) 
are  valid  only  in  the  region  which  is  0(h)  distance  away  from  the  bubble  rim. 
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In  the  thin  film  region  between  the  plates  and  the  bubble,  the  flow  field  is 
represented  by  a  Couette  flow  if  the  bubble  surface  is  rigid.  In  this  case  the  thickness  of 
the  wetting  film  can  be  given  as10-20-23-27 

f  =  1.337  -(2CaB)f     where    Ca„=^  3.7 

2  a 

Here  Can  is  the  capillary  number  based  on  the  component  of  the  bubble  velocity  normal 
to  the  bubble  surface  in  the  xy-plane,  and  a  the  equilibrium  interfacial  tension.  In  the 
absence  of  surfactant,  the  bubble  surface  is  stress-free  and  the  film  thickness  can  be 
expressed  by  the  same  equation  but  with  the  constant  2  in  the  parenthesis  replaced  by  1 . 

Equations  (3.1),  (3.5),  (3.6)  and  the  Couette  flow  along  with  the  film  thickness 
described  by  Equation  (3.7)  represent  the  three-dimensional  velocity  field  around  the 
translating  bubble  excluding  the  small  region  of  0(h)  in  the  immediate  vicinity  of  the 
bubble  rim.  Using  this  velocity  field,  an  analytic  description  of  the  translational  velocity 
of  the  elliptic  bubble  can  be  obtained  from  the  x-directional  momentum  balance  on  a 
control  volume  surrounding  the  bubble. 

An  integral  form  for  the  x-directional  momentum  balance  for  a  domain  D  which 
encapsulates  the  bubble  and  moves  with  the  same  velocity  U  as  the  bubble  may  be 
written  as  (Figure  3.1) 

ji\pvxdV  =  -\s{p[niv-\j)]vx+(np-n-T)x}dS  3.8 

This  integral  is  zero  at  steady  state  and  the  surface  integral  can  be  evaluated  using  the 


42 


V 


V 


film  region 


edge 


Bottom 


Figure  3.1       Plan  view  and  side  view  of  an  elliptic  bubble  and  a  control  volume  moving 
with  the  bubble 
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prescribed  velocity  field  over  the  four  separate  surfaces  as  shown  in  Figure  3.1;  top, 
bottom,  edge  of  the  control  volume  D,  and  the  bubble  surface. 

On  the  bubble  surface,  v  =  U  and  the  integral  over  the  bubble  surface  represents  the 
drag  force  (FD)  acting  on  the  bubble: 

|    (np-n-x)  dS  =  -FD  39 


On  the  top  and  bottom  surfaces  of  D,  the  first  term  of  the  integrand  in  Equation  (3.8) 
is  zero,  since  (v  -  U)  is  perpendicular  to  n.  The  second  term  is  evaluated  separately  for 
the  film  region  directly  above  (or  below)  the  bubble  and  the  region  which  is  0(h)  distance 
away  from  the  bubble  rim  to  give 

J       (np-n-x)  dS—^^He*  _My  3.10 

J,-*  -  r*.  )\{y  -  U)e^  +  U]\  -     W     (At) 

As  noted  in  Figure  3.1,  £=£*  and  ^0  represent  the  edge  (or  rim)  of  the  control  volume 
and  the  bubble,  respectively.  The  last  term  in  Equation  (3.10),  which  is  the  contribution 
from  the  Couette  flow  in  the  thin  film  region,  is  zero  if  the  bubble  surface  is  clean  and 
stress-free.  Here  A  is  the  plan  form  area  of  the  bubble  (i.e.,  A  =  rcab)  and  Ca  the  capillary 
number  based  on  the  bubble  velocity  U.  I  is  a  constant  given  as  a  definite  integral  which 
accounts  for  the  film  thickness  variation  in  the  transverse  direction  (i.e.,  y-direction): 
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/=jH*!-iyf[i->2]> 


3.11 


Once  k{=  a/b)  is  specified,  I  can  be  determined  numerically  and  the  results  are  given  in 
Figure  3.2  for  various  values  of  k.  When  the  bubble  is  circular  (i.e.,  k=\),I=  0.91. 

Finally  at  the  edge  of  the  control  volume,  the  only  non-zero  contribution  to  the 
surface  integral  is  due  to  the  pressure  which  is  given  as 


\(np)xdS  = 


12//  TtC1 


h     4    \k-\ 


(V-U)  +  Ve2i'  -e'U 


u+(v-u) 


k  +  l 
k-\ 


3.12 


From  Equation  (3.9),  (3.10)  and  (3.12),  the  final  form  of  the  x-directional  momentum 
balance  for  a  large  control  volume  (i.e.,  §*  -><»)  is  obtained  as 


p   =12uttc2  \     k  +  \ 
D       h     4   1     k-\ 


(V-U)  +  U 


k-\ 
k  +  l 


+ 


(V-U) 


_2 16\iU(abI)  3.13 

k~l)     2A22hcJ 


For  a  freely  suspended  bubble  with  zero  inclination  angle  (i.e.,  9  =  0  in  Fig.  1),  FD=  0  and 


Equation  (3.13)  reduces  to 


U  =  (k  +  \)V- 


2 

Ca1 


3.14 


Ca*  +0.2*7 


When  the  bubble  is  driven  by  gravity  only,  V  =  0  and  Fj)  =  -Apg(7rab/z)sin9,  since  the 
drag  force  is  balanced  by  the  buoyancy  force.  Thus, 


U  =  kU' 


2 

Ca3 


Ca*  +0.2*7 


♦      Aph  gsinO 

where     U    = 

12u 


3.15 
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Figure  3.2  Numerical  value  of  the  definite  integral,  Equation  3.11  for  various  values  of 
the  shape  factor  k. 


46 


When  the  bubble  is  driven  by  both  pressure  gradient  and  buoyancy,  Equation  (3.14)  and 
(3.15)  are  combined  to  give  a  general  expression  for  an  elliptic  bubble  as 


U  =  fi[(k  +  l)V  +  kU']        where    J3  = 


Ca 


\  3.16 


2 

Ca>+02kl 


As  Figure  3.2  indicates,  /  is  0(1)  unless  k  is  extremely  large.  Thus  (3  is  0(Co2/3). 
This  order  of  magnitude  decrease  in  the  bubble  velocity  is  due  to  the  large  drag  in  the  thin 
film  region  where  the  bubble  surface  is  rigid.  Since  this  drag  is  proportional  to  the  film 
thickness,  the  bubble  velocity  is  also  proportional  to  Co2'3.  In  the  absence  of  Marangoni 
effect,  on  the  other  hand,  the  bubble  surface  in  the  thin  film  region  is  stress  free.  In  this 
case,  the  expression  for  the  bubble  velocity  is  equivalent  to  setting  p  =  1  in  Equation 
(3.16)  since  the  last  term  in  Equation  (3.13)  is  absent.  Thus  the  result  of  Taylor  and 
Saffman  is  recovered  in  which  U  =  2V for  a  circular  bubble,  U  >  2V  and  U  <  IV  for  an 
elliptical  bubble  elongated  in  the  flow  direction  (i.e.,  k  >  1)  and  in  the  transverse 
direction  (k  <  1),  respectively. 

3.3  Comparison  with  Experimental  Results 

Retarded  motion  of  aqueous  drops  containing  sodium  dodecyl  sulfate  (SDS)  in  an  oil- 
filled  Hele-Shaw  cell  was  reported  in  the  previous  chapter  for  a  pressure  driven  flow.  In 
Figure  3.3  the  present  predictions  for  ellipses  with  various  aspect  ratios  are  compared 
with  experimental  results  for  a  drop  which  is  13  mm  in  diameter  (b/L  =  0.04).  The  drop 
contained  SDS  at  5,  10  and  20%  of  the  critical  micelle  concentration.  At  these 
concentrations  of  SDS,  the  data  sets  are  indistinguishable  from  one  another  and  show 
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reasonable  agreement  with  current  predictions  for  ellipses  that  are  elongated  in  the 
transverse  direction  (i.e.,  k  <  1).  While  the  drops  in  the  experiment  were  not  exactly 
elliptical  (Figure  2.2),  they  were  nearly  elliptical  and  elongated  in  the  transverse 
direction. 

The  present  analysis  assumes  the  bubble  shape  a  priori  and  consequently,  the 
evolution  of  bubble  shape  with  increasing  capillary  number  cannot  be  predicted. 
However,  an  interesting  observation  may  be  made  in  Figure  3.3  When  the  drop  was 
stagnant  at  Ca  =  0,  its  shape  was  circular  (i.e.,  k  =  1.0).  As  it  was  pushed  by  the 
surrounding  fluid,  the  k  value  decreased  initially  from  1.0  prior  to  exhibiting  an 
increasing  trend  with  increasing  Ca.  The  initial  decrease  in  the  k  value  occurred  at  a  very 
low  Ca  and  in  some  cases,  the  drops  simply  elongated  in  the  transverse  direction  without 
moving  until  the  average  velocity  (V)  of  the  surrounding  fluid  reached  a  certain  value. 
Unfortunately,  not  all  k  values  are  available  for  the  experimental  data  in  Figure  3.3,  and 
those  that  are  available  have  been  specified  in  the  figure.  Although  a  quantitative 
agreement  is  not  expected  between  the  observed  and  calculated  k  values,  it  is  interesting 
to  note  that  they  are  not  very  far  apart. 

It  should  be  pointed  out  that  the  retardation  factor  p  in  Equation  (3.13)  is  independent 
of  bubble  size  suggesting  that  the  same  level  of  retardation  may  occur  regardless  of  the 
bubble  size  as  long  as  their  shapes  are  identical.  One  of  the  simplifying  assumptions  in 
the  present  study,  however,  is  that  the  entire  bubble  surface  is  rigid.  While  this 
assumption  is  plausible  for  small  bubbles,  it  may  fail  for  large  bubbles  since  a  stress-free 
mobile  surface  region  may  emerge     in  the  front  part  of  the  thin  film  region20. 
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Figure  3.3    Comparison  of  experimental  results  and  the  theoretical  estimate  for  pressure 
driven  flow  (SDS  concentration:    5%  ;  A  10%;  x  20%  of  CMC) 
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Consequently,  big  bubbles  may  move  with  a  higher  velocity  than  the  present  estimate. 
Furthermore,  as  pointed  out  in  the  previous  chapter,  bubbles  which  are  large  enough  to 
have  a  mobile  surface  region  for  a  given  surfactant  concentration,  may  experience  a 
transition  at  a  certain  value  of  Ca  in  that  the  bubble  velocity  becomes  greater  than  the 
average  velocity  of  the  surrounding  fluid  (i.e.,  U>  V).  Such  a  transition,  then,  will  make 
it  possible  for  the  bubble  to  take  on  an  unusual  shape  such  as  long-  or  short-tails  as 
observed  previously.6  This  argument,  of  course,  is  only  qualitative  and  remains  to  be 
proven  rigorously. 

Experiments  with  air  bubbles  in  an  oil-filled  Hele-Shaw  cell  have  shown  that  the 
bubbles  move  with  a  velocity  comparable  to  the  predictions  of  Taylor  and  Saffman.  In 
addition,  they  are  always  elongated  in  the  flow  direction  unlike  the  bubbles  influenced  by 
surfactants.  This  difference  in  the  steady  shapes  between  clean  bubbles  and  surfactant- 
laden  bubbles  may  also  be  explained  by  the  current  study  at  least  qualitatively.  When 
two  identical  elliptic  bubbles  move  with  the  same  velocity  but  with  different  orientation 
(i.e.,  one  moving  with  its  longer  axis  in  the  flow  direction  and  the  other  with  its  shorter 
axis  in  the  flow  direction),  our  calculation  indicates  that  the  prolate  bubble  experiences 
smaller  drag  force  (or  smaller  energy  dissipation)  than  the  oblate  one  in  the  absence  of 
surfactant  effect.  In  case  of  surfactant-laden  bubbles,  on  the  other  hand,  the  prolate 
bubble  results  in  a  larger  energy  dissipation  than  the  oblate  one  due  to  the  larger  drag  in 
the  thin  film  region.  Thus,  the  elongation  in  the  flow  direction  (i.e.,  prolate  bubble)  may 
be  a  preferable  shape  for  clean  bubbles  whereas  the  oblate  shape  is  preferable  for 
surfactant-laden  ones. 
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3.4  Summary 

A  theoretical  calculation  has  been  presented  to  estimate  the  translational  velocity  of 
bubbles  or  drops  in  a  Hele-Shaw  cell  under  the  influence  of  surface  active  substances. 
Assuming  that  the  solid  wall  is  wet  by  the  surrounding  fluid  and  that  the  bubble  shape  is 
elliptical,  an  analytic  expression  has  been  derived  for  the  bubble  velocity.  The  calculated 
velocity  is  apparently  in  reasonable  agreement  with  our  own  experimental  findings  and 
also  with  those  of  Kopf-Sill  and  Homsy.  This  result  along  with  our  experimental 
findings  suggest  that  most  of  the  observations  by  Kopf-Sill  and  Homsy  are  probably  due 
to  the  influence  of  surface  active  substance  which  may  be  present  in  the  system  as 
contaminants. 


CHAPTER  4 
DISPLACEMENT  OF  A  POWER-LAW  FLUID  BY  INVISCID  FLUID 


4.1  Background 

Two  phase  displacement  flow  is  a  problem  of  industrial  importance  in  areas  such 
as  enhanced  oil  recovery,  separations  etc.,  and  has  been  studied  by  various  researchers  in 
the  past14.  When  the  displacing  fluid  is  of  low  viscosity,  then  instabilities  can  result 
causing  viscous  fingering.  Park  and  Homsy  have  theoretically  analyzed  the  two-phase 
displacement  of  Newtonian  fluids  in  a  Hele-Shaw  geometry  using  the  technique  of 
matched  asymptotic  expansions.15  Analytic  expressions  for  both  film  thickness  and 
pressure  drop  and  their  dependence  on  capillary  number  were  derived.  When  the 
displaced  fluid  is  a  shear  thinning  fluid,  these  results  can  be  very  different.  In  this  chapter, 
the  displacement  of  a  power-law  fluid  by  an  inviscid  fluid  (i.e.  air/gas)  in  a  Hele-Shaw 
cell  is  analyzed  using  the  same  technique  of  matched  asymptotic  expansions.  The  present 
analysis  is  aimed  at  deriving  the  film  thickness  and  understanding  the  stability  issue 
involved  in  this  process.  This  problem  has  practical  applications  in  the  general  area  of 
polymer  processing  and  specifically  in  applications  such  as  gas  assisted  injection 
molding. 

Figure  4.  1  is  a  schematic  of  an  inviscid  fluid  displacing  a  power-law  fluid  in  a 
Hele-Shaw  geometry.  It  is  assumed  that  the  displaced  fluid  wets  the  wall  and  there  will  be 
a  thin  film  left  behind  after  the  displacement.  The  displacement  is  slow  enough  to  ignore 
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x  =  f(z) 


Region  I       Region  II 


F(x,y,z,t)  =  0  Region  l 


>: 


2 


Fluid  2  (air/gas)  )    f  2b  Fluid  1  (Power-law  fluid) 


Figure  4.  1    Schematic  of  an  inviscid  fluid  displacing  a  power-law  fluid  in  a  Hele-Shaw 
cell 
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the  effects  of  inertia  and  the  gap  thickness  is  small  enough  to  ignore  the  effects  of  gravity. 
Under  these  assumptions,  the  interface  remains  symmetrical  about  the  center  plane.  The 
tip  of  the  interface  is  assumed  to  be  moving  with  a  constant  velocity  U  and  the  analysis  is 
carried  out  in  a  frame  of  reference  moving  with  the  tip  velocity. 

4.2  Basic  Equations 

The  governing  equations  relevant  to  the  problem  are  as  given  below. 

V.v=0  41 

Vp  =  ri(S)V2v  +  2D  •  Vrj  4.2 

Where      D  =  -(Vv  +  Vvr)  4.3a 

S  =  2(D:D)  4.3b 


n«v  =  0 
tj»T«n  =  0 

t2»T«n  =  0 
n«T»n  =  cr(V«n) 


4.3c 


n-\ 

n  =  KS  2 
K  is  called  the  consistency  factor  and  n  is  called  the  power  law  index. 

Boundary  conditions  are 

u  =  -U,    v,w  =  0  at  y  =  ±b  4.4a 


at  y  =  ±h(x,z)  44b 


Here  r|  is  the  viscosity  of  the  power-law  fluid  which  depends  on  the  shear  rate  and 
(u,v,w)  are  the  components  of  velocities  in  the  (x,y,z)  directions  respectively,  n,  t„  and  t2 


54 

are  the  unit  normal  and  unit  tangential  vectors  to  the  interface,  a  is  the  interfacial  tension 
and  T  is  the  stress  tensor.  The  origin  for  the  reference  frame  is  placed  at  the  tip  which  is 
moving  with  a  velocity  U.  b  and  h(x,z)  are  the  half  gap  thickness  and  the  location  of  the 
interface  respectively.  f(z)  is  the  projection  of  the  tip  of  the  interface  onto  the  xz  plane. 
Since  the  interface  is  symmetrical  about  the  center  plane,  only  the  bottom  half  (i.e. 
v  =  -h(x,z))is  considered  for  the  analysis.  The  following  two  dimensionless  parameters 
appear  during  the  scaling  of  the  above  problem. 

b  _       KU"b1-" 

e  =  —     and    Ca  = 

L  a 

e  is  a  parameter  which  represents  the  ratio  of  the  two  characteristic  length  scales  of  Hele- 
Shaw  geometry  and  capillary  number,  Ca,  represents  the  ratio  of  viscous  to  interfacial 
tension  forces.  A  complete  solution  to  this  problem  is  very  complicated  to  obtain,  and 
only  an  asymptotic  solution  using  the  above  mentioned  two  small  parameters  is  attempted 
in  the  following  sections. 

4.3  Scalings  and  Regions 

As  shown  in  Figure  4.  1,  the  fluid  domain  is  divided  into  three  different  regions; 

Constant  film  thickness  region  (Region  I),  front  meniscus  region  (Region  II),  and  power- 
law  fluid  region  (Region  III).  The  flow  profile  in  Region  III  can  be  obtained  by  straight- 
forward integration  of  the  governing  equations  and  is  given  in  the  next  section.  In  the 
constant  film  region  the  flow  field  resembles  plug  flow  in  the  moving  reference  frame. 
The  flow  profile  in  Region  II  can  be  understood  only  after  a  rigorous  analysis  as 
described  below.  The  following  scalings  are  used  to  non-dimensionalize  the  equations. 
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(u,v,w)~(U,U,eU) 
(x,y,z)~(b,b,L) 

p~±,    S~(U/bf 

Using  the  above  scalings,  the  Equations  4.1-  4.4  are  written  as  given  below. 


ux  +  vy  +  8  wz  =  0 


4.5 


px=Ca 


2  [      +e2(»z+wx)Sz 


4.6a 


/>v=Ca 


S2  (v„+vw+e2vs  +  — -5  2 


+e 


4.6b 


pz=Ca 


n-\ 


w„+wyy+e  w    I- 


n-\ 


n-J, 

s  2 


(uz  +  wx)Sx  + 

(vz+wy)Sy+2e2wzSz 


4.6c 


Where  5  =  l(u2x  +  v2y+  zAw]\  +  {uy  +  v^)   +  e2(uz  +  wxf  +  z2(vz  +  w^ 


at     v  =  -1, 


u  =  —  1,        v  =  vv  =  0 


4.7 


at     v  =  -/*(*,  z), 


w/7^  +  v  +  e  whz  =  0 


4.8a 


l+(V*-*.v,) 


=  0 


4.8b 
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S2 


(*i+A2) 


-  s2h:\y,  +wv) 


+  2 


=  0 


4.8c 


Ap 


_  -  ha(l  +  g2ft2)  -  g2/»g(l  +  ft,2)  +  2^*,*,  {      2Ca  S^A 


(l +  «+*>*) 


(l  +  Z^+f2/*2)  4.8d 


where 


^  =  ["  *,(«,  +  v,) -  f 2A,(v,  +  w,) -  e2hxhz(uz  +  wx)-  vy  -  h]ux  -  f4//>s] 

The  subscripts  x,y,z  in  the  above  equation  represent  partial  differentiation  of  the 
respective  variables.  To  solve  the  above  set  of  equations,  we  can  begin  by  setting  both  s 
and  Ca  equal  to  zero.  The  resulting  equations  indicate  that  pressure  is  hydrostatic  and  is 
independent  of  (x,y,z).  The  shape  of  the  interface  is  a  portion  of  a  circle  which  can  be 
obtained  by  solving  the  normal  stress  condition.  Hence  this  sub-region  (Outer  region  of 
Region  II)  may  be  referred  to  as  capillary  statics  region  where  pressure  and  interfacial 
tension  are  both  important.  In  order  to  match  the  solution  smoothly  to  the  constant  film 
region  (Region  I),  there  should  be  a  second  sub-region  (Inner  region  of  Region  II)  where 
the  interface  profile  is  smoothly  matched  to  a  constant  value.  This  may  be  called  as 
transition  region.  In  the  transition  region  both  pressure  and  viscous  effects  are  important 
and  the  variables  need  to  be  rescaled  accordingly.  Under  the  lubrication  approximation, 
the  scalings  relevant  for  this  region  can  be  derived  by  balancing  the  pressure  and  viscous 
terms  in  the  momentum  equation  and  by  balancing  the  pressure  and  interfacial  tension 
forces  in  the  normal  stress  condition.  These  scalings  are  different  from  the  ones  used  for 
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I  II  III 

A  I  B 


Figure  4.  2  Schematic  showing  the  scales  and  coordinate  systems  for  different  regions. 
(A:  Transition  region;  B:  Capillary-statics  region) 
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the  Newtonian  case  and  they  depend  upon  the  power-law  index  and  capillary  number  as 
mentioned  below. 


(i^-(.*3T.» 


\xtytz)  = 


CaK 

x  +  l    y  +  \ 

~cJ,Ca2k' 


where  k  = 


1 


2n  +  \ 


h   =  -L±;     p  =  p;     S  =  Ca4kS 
La 


The  variables  in  the  transition  region  are  denoted  by  an  overbar.  As  indicated  in  Figure 
4.2,  the  origin  for  the  transition  region  is  x  =  -I  and  v  =  -1.  The  unknown  /  will  be 
determined  later  using  the  matching  conditions.  Using  the  newly  scaled  variable,  the 
governing  equations  valid  for  the  transition  region  are  given  below. 


ux  +vy  +e2Cakwz  =  0 


4.9 


P~x       = 


S  2 


rCa2kUxx 
[+e2Ca4kUv 


n-3 


+  a^lsi 


2ca2kUxSx 


(uy  +  Ca2k Vx\Sy 
e2Cau{Caku-z+Wx)s-Z/ 


4.10a 


—  _n-i 

p-  =  Ca2kS2 


'Ca2kv^ 


+Vyy 

+e2Ca4kvl 


n-3 


+  n^Ca2kS~2 


i^-y  +  Ca2kVx^jSx 

+    2VyS~y 

+e2CakS-2(Ca3kv-:+Wy)) 


4.10b 


n-l 


'Ca2kwlx 


Cakp-=S  2 


+  Wyy 

+z2Ca4kwTz 


n-3 


+      ^S2 


' Ca2k(Cak uz  +  wx)s-x  + 
(Ca3kv'z  +Wy)~S^y 


\         4.10c 


+2z2Ca4kw-zS: 
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Where  S  = 


2Ca2k(ul  +  vly+  e4Ca2k  wfj  +  [uy  +  Ca2k  v~x) 
+z2Ca2k(CakU:  +w~xf  +s2(Ca3*vi  +  w~yf 


at    y  =  0, 


u-  -1,        v  =  w  =  0 


4.11 


at    y  =  h(x,z), 


2^„k~ 


uhx  -  v  +  s  Ca  wh~z  =  0 


4.12a 


n-l 

S2  <! 


+Ca2*(Ca*«;+^)(S-E2Ca2*/jf)     +     Ca3*(fcv;  -fcvi) 


=  0 


4.12b 


n-l 

5  2 


V      X  Zy 


hx[uy+Ca2kvx 


+  E2Cakh~z\Ca3kv-z  + 


) 


■  lca3kvz+wy\ 


+  2Ca 


2k 


h2-(uX-V-y) 


+  s2Ca2k  h-U2Cak  wz  -v. 


3~kT- 


+  e  CaKhxhz\Ca    uz  +  wx 


k  - 


-hxiuy+Ca2kvx\-s2Cakh~zICa3kvz+Wy\ 


=  0 


4.12c 


^p  = 


\  +  £2Ca4kh2J+€2Ca2kh7z[\  +  Ca2khlj-2e2Ca4khxhzhx-z 


+ 


2Ca-2kn+k+\  s   2 


\  +  Ca2kh2-+e2CaUh2^2 


n-l 


M2d 


{\  +  Ca2kh2-  +  e2CaUh?) 


hx(uy+Ca2kvx)  +  £2Cakhz(Ca3kvz~  +  wy^\ 


-e2Ca3khxhz(cakuz+wx 
-Ca2kh-ux  -  £4Ca5kh-w-z 

X  z 


j-Vy 
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4.4  Small  Parameter  Expansions 

To  obtain  the  solution,  all  the  variables  (u,  w,  p,  h  etc..)  are  expanded  as  given 

below. 

p  =  (p00  +Caap01  +Ca2ap02  +..)  +  e2(p10  +  CaapU  +Ca2apn  +...)  +  ...     4.13 

The  choice  of  gauge  function,  a  ,  is  obtained  by  comparing  all  the  exponents  that  appear 
in  the  governing  equations  and  boundary  conditions  and  choosing  the  least  common 
multiple  of  them  so  that  matching  is  possible  at  every  order.  For  example,  the  choice  of  a 
few  gauge  functions  are  as  given  below. 

a    =    Ink    if  «  =  -;-;  -etc.  414 

2  4    6 

,      >r        111 
=     nk     it  »  =  —;—;—  etc. 

3  5    7 

Since  the  gauge  function  depends  on  the  power-law  index,  only  the  leading  order  solution 

is  attempted  in  the  following  sub-sections.  Once  the  choice  of  gauge  function  is  known, 

higher  order  corrections  can  be  obtained  following  the  procedure  of  Park  and  Homsy. 

Consequently,  in  the  present  analysis  the  dynamic  pressure  jump  across  the  interface  can 

not  be  determined  since  it  appears  at  a  higher  order. 

4.4.1     Capillary  Statics  Region  : 

Substituting  4.13   into  4.5-4.8,   the   leading   order  governing   equations   and 

boundary  conditions  for  the  capillary  statics  region  are  as  given  below. 

„°°  _■_  „00     n 

ux  +vy  =0  4_13 


„00  _      00  00       ft 

Px    -Py    =Pz    =0  4.14 
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Boundary  conditions,      w00  =  -1,        v00  =  w°°  =  0  at  y  =  -1  4-15a 

at     y  =  -h00(x,z)      u00h°x0  +  v00  =  0  4-15b 

An00  =       ~h™  4,15c 


(l^f) 


[(W  -h?w?)  +  2h?h?u?  +hf(u°°  +w?)+(h?v?  -^00vz00)]  =  0  4.15d 

\-h?(u?  +  vf  )  +  2/,f>»  -  v»)  + Af  («»  +v»)|  =  0 

4.15e 

As  evident  from  the  above  equations,  the  viscous  force  term  in  the  momentum  equations 
is  negligible  and  the  pressure  is  independent  of  (x,  v,  z).  The  shape  of  the  interface  can  be 
obtained  by  integrating  the  normal  stress  balance  and  is  given  as  below. 


Ap 


oo 


1 


w/wH'  *  416 


The  two  boundary  conditions  that  were  used  in  the  integration  and  which  apply  at  all 
orders  are 

hx-+-co        as      x->/(z)  417 

h  =  0  at      x  =  f(z) 

The  solution  4.16  indicates  that  the  tip  of  the  interface  is  a  portion  of  a  circle  modulated 

in  the  z-direction  by  f(z).  This  solution  contains  one  unknown  Ap00,  which  will  be 

determined  using  the  matching  conditions.    Equation  4.16  is  not  uniformly  valid  in  the 

entire  region  and  it  needs  to  matched  with  the  constant  film  region  (Region  I). 
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4.4.2     Transition  Region 

Variables  in  Equations  4.9-4.12  are  expanded  using  4.13  and  the  leading  order 
terms  are  gathered  to  obtain  the  governing  equations  for  the  transition  region. 

4.18 


-00      -00       n 
Ux    +  V  v    =0 


oo      I  -o?" 


pT  =  o 


-00 
A        ~ 59-  /  ,\  Uyy  —00 

0=Wyy        +        (n-\)-=WWy 

Uy 


at   v  =  0 


oo, 


at  v  =  h    (x,  z) 


-oo       ,         -oo     —oo     n 
u     =-1,         v     =  w     =0 

-ooToo     -oo     n 
u    hx  -  v     =0 

—00       Too 
Ap     =  hxx 


uy  =0 


4.19a 
4.19b 

4.19c 

4.20 
4.21a 
4.21b 
4.21c 


By  solving  these  equations  and  matching  with  the  constant  film  thickness  region  (Region 
I)  the  velocity  profile  can  be  determined.  Once  the  velocity  field  is  known,  then  normal 
stress  along  with  kinematic  condition  results  in  an  evolution  equation  for  the  shape  of  the 
interface  as  given  below. 


t0Q__ 

flxxx  — 


2n  +  \ 


Too     -00 
h     -t 


h 


2»+l 

00~"~ 


4.22 


When  n  is  set  equal  to  one,  the  above  equation  reduces  to  the  Newtonian  case  (Park  and 


-00 


Homsy  1984).   t      represents  the  leading  order  approximation  for  the  constant  film 
thickness  which  will  be  determined  later  using  the  matching  conditions.   For  numerical 
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integration,  the  above  equation  is  transformed  into  a  canonical  form  by  using  the 
following  transformations. 


-oo 
H00  =  —  4.23 

-00 


x+s  n+2 

X  =  TT^B      where  P  =  — r~  4.24 


s  is  an  arbitrary  shift  factor  which  is  determined  by  higher  order  matching  conditions. 
With  these  transformations,  the  evolution  equation  is  written  as  follows. 


»(H00-\)" 


oo       (2n  +  \Y\"     -V  4.25 

Hxxx~{     n    )    lHw\^ 

The  condition  for  matching  with  region  I  now  becomes 

//oo  _>,         ^      ^^.^  4.26 

Equation  4.25  can  be  numerically  integrated,  using  small  slope  and  small  curvature  as 
initial  conditions.  A  4,h-order  Runge-Kutta  routine  is  used  to  do  the  numerical  integration. 
The  shape  of  the  interface  depends  on  the  power-law  index  and  it  is  plotted  in  Figure  4.3 
for  the  case  where  n  =  0.5  along  with  the  profile  for  a  Newtonian  case  (i.e.  n  =  1.0).  Since 
as  X  ->  +oo  the  interface  profile  H00  also  goes  to  infinity,  it  is  possible  to  approximate  the 
profile  using  a  quadratic  form  as  given  below. 

//00(x,z)    =     l-CX+C,(z)X  +  C2(z)  4-2? 

Equation  4.27  can  be  written  in  terms  of  the  original  scalings  as  given  below. 
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Figure  4. 3   Comparison  of  interface  profiles  for  Newtonian  («=1)  and  power-law 
(n=0.5)  fluids. 
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-oo 


h    (x,  z)     = 


X    + 


oo2A-i 


It 


C..s 


-oo2^-'       -oo"-1 
^2/  t 


cy     c.s    _  V 

V2r  t 


4.28 


-uo 


The  unknown  /    will  be  determined  by  the  following  matching  conditions. 

4.4.3  Matching  Conditions 

The  solution  in  the  transition  region  has  to  be  matched  with  the  solution  in  the 
capillary  statics  region  using  the  following  equation. 

4.29 


Urn  \\-Ca2kh(x,z)}=   Hm  h(x,z) 


This  limits  are  to  be  interpreted  in  terms  of  the  matching  principle  by  Van  Dyke  (1964). 
By  expanding  h(x,z)  about  x  =  -/  using  Taylor  series  expansion,  rewriting  the  expansion 
in  inner  variables  and  comparing  it  term  by  term  with  the  left-hand  side,  matching 
conditions  for  each  order  can  be  determined  as  given  below.  Higher  order  matching 
conditions  can  be  derived  once  the  choice  of  a  is  known. 


0{e°Ca°) 
0(s°Cak) 
0(z°Ca2k) 


h°°(-l°,z)  =  \ 
A?(-/V)  =  0 

«(-/V)=-ce(?°°) 


4.30 
4.31 


2w+l 
3 


4.32 


Physically  these  conditions  mean  that  the  outer  static  solution  meets  the  wall  with  zero 
slope  and  the  curvature  matches  with  that  of  transition  region  solution.  Solving  the  above 
equations,  the  following  unknowns  in  the  solution  can  be  determined. 
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A/0  =  1 


/°=l-/(z) 


-00 


3 
2n+l 


4.33a 
4.33b 
4.33c 


The  above  solution  determines  the  location  of  the  origin  and  also  the  film  thickness  to  a 

leading  order.    Since  the  film  thickness  is  stretched  by  Ca2k ,  the  dimensionless  film 
thickness  is  given  by 


,°°=(C0)^ 


KU"b 


hlI-a 


2w+l 


4.34 


Figure  4.4  displays  the  dependence  of  the  constant  Co  on  the  power-law  index  n.  As 
n  -» 1  the  constant  C0  approaches  1.337  which  is  the  value  for  the  Newtonian  case.  In 
Equation  (4.34),  if  K  is  replaced  by  u0/y0  where  n0is  apparent  viscosity  and  y0  is  a 
reference  shear  rate,  then  the  film  thickness  is  given  as  follows. 


(oo 


(c.) 


3_ 
2n+\ 


fr/oHo 


.2(H 


In 


4.35 


'Ca'2-1     where  Ca'=^- 


Figure  4.5  is  a  plot  of  Equation  (4.35)  which  is  a  straight  line  on  a  log-log  scale.  All  the 
variables  influence  the  intercept  while  only  the  power  law-index  influences  the  slope  of 
the  straight  line  plot.  In  the  Newtonian  case,  all  the  fluids  have  a  master  curve  (i.e.  the 
plot  for  n  =  1)  whereas  each  power-law  fluid  has  a  separate  curve  uniquely  defined  for  the 
given  operating  conditions. 
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4.5  Stability  Analysis 

In  order  to  carry  out  the  stability  analysis,  it  is  important  to  derive  the  flow  field  in 
the  power-law  fluid  region.  Although  far  away  from  the  interface,  the  flow  is  only  along 
the  x-direction,  near  the  interface  a  two  dimensional  flow  field  is  considered.  The  solution 
to  this  region  is  fairly  straight  forward.  It  is  assumed  that  the  gradient  of  components  of 
velocities  in  x,z  direction  are  negligible  and  only  the  gradients  in  the  gap  direction  (i.e.  v 
direction)  are  important.  The  simplified  governing  equations  and  boundary  conditions  are 
given  below. 


PX  =  K 


n-l 
S  2  1 


where    S  —  uy  +  wy 
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4.37a 
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The  above  set  of  equations  can  be  solved  to  obtain  the  velocity  field  as  given  below. 
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The  components  of  velocity  are  averaged  in  the  gap  direction  (i.e.  y  direction)  and  the 
resulting  velocity  profile  is  written  in  the  vector  form  as  given  below. 


v"*  =  -■ 


nb 


2«  +  l 


(b\n       lr7,i^„ 

UJ       N  ■  *P  4.40 


When  n  =  1,  the  above  equation  reduces  to  Darcy's  law.  However,  in  the  case  of  power- 
law  fluids,  it  may  be  noted  that  the  average  velocity  field  is  non-Darcian.  The  interface  is 
assumed  to  be  flat  across  the  gap  and  moving  with  a  constant  velocity  U.  Then,  the 
average  velocity  in  the  power-law  fluid,  uavg,  can  be  related  to  the  velocity  of  the 
interface,  U,  as  given  below.  The  average  velocity  in  the  lateral  direction  is  zero. 


ave     tt~          i         ~     b  -  film  thickness  4.41 

if9*  m  U  t        where   t  = 


The  location  of  the  interface  is  denoted  by  x  =  r\0  which  is  a  function  of  time  and  without 
loss  of  generality  the  pressure  inside  the  inviscid  fluid  may  be  set  equal  to  zero.  Now, 
equation  (4.41)  is  substituted  in  equation  (4.40)  and  the  steady  state  solution  to  the  depth 
averaged  pressure  field  can  be  obtained  as  given  by  the  following  equation. 

s_    ((2n  +  l)7uV 


{         nb 


K(  \ 

J{X-^  4.42 


Equations  (4.41)  and  (4.42)  constitute  the  base  state  solution  or  the  steady  state  solution 
for  the  stability  analysis.  The  stability  of  this  solution  is  analyzed  following  the  standard 
procedure  for  the  linear  stability  analysis.  The  steady  state  velocity  and  pressure  field  is 
perturbed  by  introducing  disturbances  in  all  the  physical  variables  as  given  below. 
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u  =  U  +  u' 


W=W'  P  =  PS+P'  X  =  TJ0  +  Tj'(z,t)  4 


43 


Here  the  superscript  prime  indicates  that  it  is  a  perturbation  from  the  steady  state  solution. 
By  substituting  in  equation  (4.40),  the  relation  between  the  perturbed  velocity  and 
pressure  can  be  written  as  follows. 


!/'=- 


w  =  - 


4.44a 


4.44b 


According  to  the  linear  stability  theory,  the  perturbation  in  the  velocity  field  can  be 
written  as  follows. 


u'  =  V|/(.x:)     exp  \ikz  +  8/1 


4.45 


Here  k  is  the  wave  number  and  8  is  the  growth  rate  of  the  disturbance.  Using  equation 
(4.45)  and  the  continuty  equation  (i.e.  u'x  +  w'z  =  0),  an  expression  for  w'  can  be 
obtained.  By  substituting  u',  w'  in  equation  (4.44)  and  eliminating  pressure  the  following 


ordinary  differential  equation  is  obtained. 


\\i"  -nk  \\i  =  0 


4.46 


Far  away  from  the  interface,  into  the  power-law  fluid  region,  the  disturbance  dies  out. 
Hence  the  above  differential  equations  can  be  solved  assuming  that  as  x  — »  oo,  *P  — »  0 . 
The  resulting  solution  is  as  given  below. 


i|/(jc)  =  a  expl  -ynk2  x 


4.47 


Now,  the  components  of  the  perturbed  velocity,  u',w' ,  are  determined  to  a  constant.  By 
substituting  them  in  equation  (4.44)  and  integrating,  the  pressure  field  in  the  power-law 
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fluid  region  can  be  determined.  As  indicated  in  equation  (4.43),  the  instantaneous 
location  of  the  interface  is  given  by  x,  which  can  be  differentiated  to  give  instantaneous 
velocity  in  the  x-direction  as  mentioned  below. 

(x  -  t]0)t  =  n't  =  u'\        ^  ,  4.48 

Equation  (4.48)  can  be  integrated  with  respect  to  time  to  obtain  an  expression  for  the 

perturbed  location  of  the  interface  (i.e.  rj' ).  The  pressure  jump  across  the  interface  can  be 

written  as 

,       ,      n       ,  4.49 

P\-P2=-OVizz 

The  constant  n/A  appearing  in  equation  (4.48)  was  first  derived  by  Park  and  Homsy 
(1984)  as  an  0(s2)  correction  that  accounts  for  the  curvature  change  in  the  z-direction. 
Substituting  the  deviation  variables  in  the  above  boundary  condition,  the  dispersion 
relation  for  the  growth  rate  of  the  disturbance  can  be  written  as  follows. 


.     Uk 


.  _  J         nb        I     ban    2 
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4.50 


When  n  is  set  equal  to  one,  the  dispersion  relation  for  the  Newtonian  case  is  recovered 
(Saffman  &  Taylor  1958).  Figure  4.6  is  a  graph  of  5  as  a  function  of  the  wave  number.  If 
8  is  positive  then  the  disturbance  grows  at  an  exponential  rate  and  the  interface  is 
unstable  to  small  disturbances.  If  8  is  negative,  then  the  disturbance  is  damped  at  an 
exponential  rate  and  the  motion  is  stable  to  small  disturbances.  As  is  evident  from  the 
graph,  there  is  an  optimum  value  for  the  wave  number  for  which  the  growth  rate  is 
maximum  and  beyond  which  value,  the  growth  rate  monotonically  decreases.  So  8  =  0 
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represents  the  cut-off  wave  number  whose  expression  is  as  given  below. 


\AK  I(2n  +  1)TU]  4.51 

cut-°ff=pob\        nb       J 

Hence  the  wave  numbers  greater  than  the  cut-off  wave  numbers  die  out  while  the  smaller 
ones  grow  at  different  growth  rates.  The  wave  number  with  the  maximum  growth  rate  can 
be  obtained  by  differentiating  8  with  respect  to  k  and  setting  it  equal  to  zero  and  it  is 
related  to  the  cut-off  wave  number  as  given  below. 

km=^f-  4.52 

In  equation  (4.51),  the  consistency  index  K  can  be  replaced  by  apparent  viscosity  (u0)  and 

a  reference  shear  rate(y0)  as  K  =  ju0/7o~l  >  men  equation  (4.52)  can  be  plotted  for 
various  power-law  fluids  as  given  in  Figure  4.7.  While  making  this  comparison,  it  should 
be  kept  in  mind  that  the  reference  shear  rate  should  be  high  enough  so  that  power-law 
model  is  applicable  to  all  the  fluids  being  compared.  Because,  one  of  the  limitations  of 
power-law  model  is  that  it  can  not  describe  the  viscosity  behavior  at  low  shear  rates. 

The  result  in  Figure  4.7  indicates  that,  the  wave  number  for  maximum  growth  rate 
(kmax)  decreases  as  the  power-law  index  reduces.  In  other  words,  the  wave  length  for  the 
maximum  growth  (X,max)  is  greater  for  power-law  fluids.  To  experimentally  observe  this 
wave  length,  the  width  of  the  Hele-Shaw  cell  (L)  should  be  larger  than  A.max.  Hence,  for 
fluids  with  similar  viscosity  at  a  constant  reference  shear  rate,  it  may  be  concluded  that 
the  interface  will  be  less  sinuous  for  a  power-law  fluid  than  for  a  Newtonian  fluid. 
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4.6  Summary 

The  displacement  of  a  power-law  fluid  by  an  inviscid  fluid  is  analyzed  and  a 
leading  order  solution  is  obtained  using  the  technique  of  matched  asymptotic  expansions. 
For  slow  displacement,  the  shape  of  the  interface  in  the  capillary  statics  region  is  a 
portion  of  a  circle  which  is  smoothly  matched  to  a  parabolic  profile  in  the  transition 
region.  Under  the  lubrication  approximation,  the  shape  of  the  interface  in  the  transition 
region  is  obtained  by  numerically  integrating  the  evolution  equation.  Matching  conditions 
with  the  capillary  statics  region  determine  the  film  thickness  of  Region  I.  The  flow  field 
in  Region  III  is  determined  and  found  to  be  non-Darcian.  Assuming  that  the  interface  is 
flat,  the  steady  state  solution  for  the  velocity  and  pressure  field  is  determined.  Stability 
analysis  for  this  steady  state  solution  is  carried  out  to  determine  the  wave  number  with 
the  maximum  growth  rate.  The  result  indicates  that  the  interface  will  be  less  sinuous  for 
power-law  fluids  than  Newtonian  fluids. 


CHAPTER  5 
CONCLUSIONS 


The  experimental  findings  of  Kopf-Sill  and  Homsy  on  the  motion  of  air  bubbles 
in  a  Hele-Shaw  cell  filled  with  Glycerin-water  mixture  have  raised  the  following  issues. 

•  The  velocities  of  small  air  bubbles  were  found  to  be  an  order  of  magnitude  smaller 
than  the  theoretical  predictions  of  Saffman  and  Taylor. 

•  A  variety  of  unusual  steady  bubble  shapes  and  a  sudden  transition  of  bubble  velocity 
and  shape  at  a  critical  capillary  number  were  observed. 

The  present  study  explains  the  above  observations  in  the  context  of  surface  active  organic 
contamination. 

To  study  the  influence  of  surfactants,  experiments  were  conducted  both  using  air 
bubbles  and  water  drops  containing  surfactant  in  an  oil-filled  Hele-Shaw  cell.  While  the 
motion  of  air  bubbles  is  in  agreement  with  the  predictions  of  Saffman-Taylor  theory,  the 
velocity  of  water  drops  is  an  order  of  magnitude  smaller  than  the  predicted  value.  Unlike 
the  air  bubbles,  the  water  drops  exhibited  a  variety  of  shape  transitions  similar  to  the  ones 
observed  by  Kopf-Sill  and  Homsy.  In  the  presence  of  surfactants,  the  bubble  motion  can 
induce  a  surface  concentration  variation  resulting  in  Marangoni  stresses  along  the  drop 
interface.  These  stresses  can  cause  the  drop  surface  to  be  rigid.  In  the  presence  of  a 
wetting  film,  the  dissipation  in  the  thin  film  region  can  be  significantly  high  and 
consequently  the  velocity  of  the  water  drops  is  decreased  by  an  order  of  magnitude.  The 
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observed  shapes  can  also  be  explained  using  the  same  argument.  The  motion  of  the 
bubble  creates  a  spatial  variation  of  surface  tension.  In  order  to  keep  the  pressure 
constant,  the  curvature  normal  to  the  gap  direction  adjusts  itself  resulting  in  a  variety  of 
shapes.  Theoretical  modeling  was  done  for  an  elliptic  bubble  assuming  that  the  entire 
surface  is  rigid.  The  theoretical  predictions  agree  with  the  experimentally  observed 
velocities.  This  calculation  also  suggests  that  a  sideways  elongated  shape  is  the  preferred 
one  in  the  presence  of  surfactants. 

In  Chapter  4,  the  displacement  of  a  shear  thinning  power-law  fluid  by  an  inviscid 
fluid  is  analyzed.  The  flow  domain  is  divided  into  three  different  regions  and  a  leading 
order  solution  is  obtained  using  the  technique  of  matched  asymptotic  expansions.  Since 
the  displacement  is  very  slow,  in  the  capillary  statics  region  the  shape  of  the  interface  is 
nearly  circular.  Under  the  lubrication  approximation,  the  interface  shape  in  the  transition 
region  is  obtained  by  numerically  integrating  an  evolution  equation.  Matching  condition 
with  the  capillary  statics  region  provides  an  expression  for  the  film  thickness  in  the 
constant  film  thickness  region.  In  the  case  of  Newtonian  fluids,  irrespective  of  individual 
physical  properties,  capillary  number  uniquely  determines  the  film  thickness.  On  the 
other  hand,  in  the  case  of  shear  thinning  power-law  fluids,  the  relationship  is  unique  to  a 
given  fluid  and  flow  conditions.  Just  as  in  the  Newtonian  case,  as  the  Capillary  number  is 
increased,  the  film  thickness  increases  linearly  on  logarithmic  scale.  However,  at  a  given 
capillary  number,  fluids  with  a  lower  power-law  index  have  a  smaller  film  thickness.  The 
flow  field  away  from  the  interface  in  the  power-law  fluid  region  is  also  solved  assuming 
that  the  gradients  of  components  of  velocities  are  negligible  in  the  directions  normal  to 
the  gap  direction.  Unlike  the  Newtonian  case,  the  depth  averaged  velocity  field  is  found 
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be  non-Darcian.  A  stability  analysis  is  carried  out  for  a  flat  interface  and  a  dispersion 
relation  is  derived  for  the  growth  rate  as  a  function  of  wave  number.  The  result  indicates 
that  the  wave  number  for  the  maximum  growth  rate  decreases  as  the  power-law  index 
decreases.  In  other  words,  the  dominant  wave  length  increases  as  the  fluid  becomes  more 
and  more  non-Newtonian.  Hence,  the  interface  will  be  less  sinuous  for  the  power-law 
fluid 
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